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Preface 


The  Aerial  Refueling  Branch  of  the  Airlift  Manage¬ 
ment  Division,  Headquarters,  Military  Airlift  Command 
(HQ  MAC / DOOR ) ,  has,  for  some  time,  required  an  analytical 
method  for  optimizing  the  use  of  aerial  refueling  on  stra¬ 
tegic  airlift  missions.  This  study  was  conceived  to  satisfy 
this  requirement.  Up  till  now,  MAC  planners  have  dealt 
with  the  problem  of  how  to  efficiently  employ  the  use  of 
tankers  either  from  past  experience  or  by  iteratively  plan¬ 
ning  a  mission  to  test  various  scenarios.  This  study  pro¬ 
vides  a  model  that  will  determine  optimal  rendezvous  points, 
fuel  offloads,  and  tanker  departure  bases,  using  the  total 
fuel  consumed  by  both  airlifter  and  tanker  as  the  measure 
of  effectiveness.  Use *of  this  model  will  insure  optimal 
use  of  tanker  resources,  save  considerable  time  during 
aerial  refueling  mission  planning,  and  will  result  in  a 
considerable  fuel  savings  to  the  Air  Force. 

Special  thanks  goes  to  Major  Gerald  Armstrong,  my 
thesis  advisor,  for  initially  providing  a  possible  solution 
technique  for  this  very  complex  problem  and  his  continual 
guidance  during  the  preparation  of  this  report.  I  would 
also  like  to  express  my  sincere  appreciation  to  Professor 
Dan  Reynolds,  my  reader,  for  his  valuable  comments,  and 


his  wife,  Phyllis,  for  her  outstanding  work  in  preparing 
the  final  copy  of  this  report.  LtCol  Tom  Baars  and  Capt 
Mark  Kahley,  both  of  HQ  MAC/DOOR,  also  provided  invaluable 
assistance  in  insuring  that  the  aerial  refueling  optimiza¬ 
tion  model  was  a  valid  one  and  that  it  could  in  fact  be 
used  operationally.  Finally,  I  want  to  thank  my  wife, 
Lenita,  for  her  faithful  support  and  sacrifice  during  some 
busy  and  difficult  times. 
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Abstract 


With  the  advent  of  strategic  airlift  in-flight 
refueling,  the  role  of  strategic  mobility  in  supporting 
our  national  policy  has  taken  on  new  dimensions.  Airlift 
of  important  military  equipment  to  any  part  of  the  world 
is  no  longer  dependent  on  enroute  bases  for  refueling. 

Also,  our  response  time  to  any  contingency  far  from  home 
can  now  be  measured  in  terms  of  hours  instead  of  days. 

This  increased  demand  on  our  limited  aerial  refueling 
resources,  plus  the  high  cost  and  doubtful  availability 
of  jet  fuel,  require  that  any  deployment  using  aerial  refuel¬ 
ing  be  done  as  efficiently  as  possible. 

The  main  objective  of  this  study  is  to  define  and 
develop  an  operationally  usable  optimization  model  that 
will,  given  airlifter  route  of  flight,  takeoff  fuel  and 
gross  weight,  enroute  wind  factors,  tanker  departure  bases, 
and  available  enroute  diversion  bases,  determine  the  aerial 
refueling  points,  fuel  offload,  and  tanker  departure  bases 
to  minimize  total  fuel  consumed  by  both  airlifter  and 
tanker.  The  key  to  this  study  is  the  requirement  that  the 
model  be  operationally  usable.  Therefore,  the  model  is 
designed  to  allow  input  of  any  mission  scenario  using 
geographic  coordinates. 
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The  various  subobjectives  presented  in  this  paper 
are  organized  with  a  view  towards  validating  and  verifying 
the  model's  output.  However,  some  questions  are  also 
answered  regarding  the  nature  of  the  aerial  refueling 
optimization  problem.  An  analysis  of  the  experimental 
design  shows  that  it  is  extremely  difficult  if  not  impos¬ 
sible  to  develop  general  "rules  of  thumb"  to  minimize  total 
fuel  consumed,  primarily  because  of  the  complex  interaction 
between  the  control  variables  within  the  aerial  refueling 
system.  It  is  found  that  cargo  load,  and  hence  airlifter 
gross  weight,  has  more  effect  on  total  fuel  consumed  than 
takeoff  fuel.  Also,  it  is  found  that  the  objective  of 
minimizing  total  fuel  consumed  is  not  always  consistent 
with  that  of  minimizing  the  use  of  tankers. 

The  model  uses  the  dynamic  programming  technique 
during  the  solution  process.  In  addition,  the  "Complex" 
Method  of  Box  is  incorporated  to  optimize,  or  minimize, 
stage  return  functions. 

It  is  anticipated  that  the  use  of  this  model  to 
plan  operational  aerial  refueling  missions  will  result  in 
significant  fuel  savings  to  the  Air  Force. 


AN  AERIAL  REFUELING  OPTIMIZATION  MODEL 


APPLIED  TO  STRATEGIC  AIRLIFT 


I  Introduction 


Background 

A  recent  SALT  negotiator  remarked  that  the  Soviets 
seemed  envious  of  only  one  American  military  capability 
the  routine  and  seemingly  effortless  way  we  refueled, 
at  night  and  in  all  weather,  our  tactical  and  strategic 
forces.  For  some  reason,  that  capability  has  not 
received  the  same  appreciation  in  the  Pentagon's  E-Ring 
these  past  several  years  [Ref  20:87]. 

The  above  quote  from  General  T.  R.  Milton,  USAF 
(ret.),  highlights  the  importance  of  our  strategic  aerial 
refueling  capability  and  notes  that  planned  substantive 
improvements  to  our  aerial  tanker  fleet  are  not  receiving 
the  proper  emphasis.  As  we  see  tighter  fiscal  policies 
evolve  in  response  to  the  ever-mounting  national  debt,  it 
is  evident  that  our  aerial  refueling  capabilities  will  not 
keep  up  with  the  increasing  number  of  aerial  refueling- 
capable  aircraft  in  the  Air  Force  inventory.  Two  very 
important  improvements  to  our  aerial  refueling  capability 
have  been  either  cut  from  planned  levels  or  delayed--the 
new  KC-10A  and  the  KC-135A  re-engining  program  are  both 
victims  of  the  proposed  Fiscal  Year  1982  budget  (Ref  10:19) 
The  increasing  demand  for  limited  aerial  refueling 
resources,  plus  the  continuing  problem  of  jet  fuel 


availability  and  cost,  demand  that  these  resources  be  used 
as  efficiently  as  possible.  This  study  will  define  and 
develop  a  computer  model  that  can  be  used  operationally  to 
optimize  the  use  of  our  tanker  fleet.  Although  this  study 
specifically  addresses  the  use  of  tankers  by  the  Military 
Airlift  Command  (MAC) ,  other  Air  Force  users,  such  as  the 
Strategic  Air  Command  (SAC) ,  will  find  the  results  equally 
adaptable  to  their  aerial  refueling  mission. 

Although  aerial  refueling  was  first  demonstrated 
feasible  in  1929  by  the  legendary  flight  of  Ira  Eaker, 
Tooey  Spaatz,  and  others,  its  full  impact  was  not  realized 
until  the  Vietnam  War  (Ref  20:87).  In  Thailand,  which  saw 
the  daily  rendezvous  of  North  Vietnam-bound  fighters  with 
an  aerial  fleet  of  KC-135s,  we  began  to  better  understand 
and  develop  the  full  potential  of  this  capability.  Aerial 
refueling  had  been  used  for  many  years  in  a  training  role; 
however,  this  was  the  first  time  the  concept  had  been 
proven  in  actual  combat.  It  also  showed  that  aerial 
refueling  had  very  important  applications  in  other  areas 
besides  the  traditional  strategic  bomber  support  role. 

Later,  during  the  Yom  Kippur  War  of  1973,  a  new 
important  role  for  aerial  refueling  began  to  emerge.  The 
ensuing  aerial  resupply  of  Israel  by  the  United  States 
showed  that  strategic  mobility  may  very  well  be  the  key 
to  success  during  a  contingency  involving  our  interests 
far  from  home.  Operation  Nickel  Grass,  the  code  name 
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assigned  to  the  Military  Airlift  Command's  resupply  mis¬ 
sion,  was  extremely  successful,  resulting  in  the  delivery 
of  22,395  tons  of  military  equipment  in  30  days  (Ref  1:122) 
It  did  serve  to  point  out  one  very  disturbing  fact,  how¬ 
ever,  and  that  was  the  dependency  of  any  global  airlift 
contingency  on  the  availability  of  enroute  bases  for  refuel 
ing.  Due  to  the  political  implications  of  any  future  con¬ 
flict,  landing  rights  in  foreign  countries,  even  for  only 
refueling,  may  be  denied.  Specifically,  if  Lajes  AB, 
Azores,  had  not  been  available  for  our  use  during  Nickel 
Grass,  another  refueling  and  staging  base  would  not  have 
been  available.  For  this  reason,  MAC  placed  its  highest 
priority  on  developing  its  already  existent  (but  dormant) 
aerial  refueling  capability. 

Today,  MAC  has  total  strategic  airlift  aerial 
refueling  (AR)  capability.  C-5As  have  been  flying  opera- 
tional  AR  missions  since  1975.  C-141Bs  have  been  flying 

training  and  operational  AR  missions  since  the  beginning 
of  1980,  and  fleetwide  modification  of  C-14lAs  to  the  AR 
version  will  be  completed  in  1982.  The  new  C-17  will  also 
have  AR  capability.  Future  limitations  on  mobility  will 
not  be  dependent  on  having  enough  AR-capable  airlifters, 
but  rather  on  how  the  use  of  our  airlift  resources  will 
be  planned.  Competition  from  other  commands  for  the  use 
of  limited  tanker  support  will  require  these  tankers  to 
be  efficiently  allocated  and  used. 


MAC  is  presently  without  a  Command  policy  specify¬ 
ing  which  missions  should  be  aerial  refueled  during  a  con¬ 
tingency  (Ref  13) .  Consequently,  contingency  planning 
with  analytical  techniques  is  done  by  using  AR  only  on  a 
random  basis.  One  part  of  the  problem  lies  in  the  fact 
that  MAC  planners  do  not  have  a  means  of  determining  the 
optimal  use  of  tankers.  Another,  and  more  basic,  problem 
is  that  SAC  has  not  made  available  to  MAC  planners  the 
number  of  tankers  that  will  be  available  to  them  in  the 
event  of  a  contingency  (Ref  13) .  An  optimizing  AR  plan¬ 
ning  tool  is  first  needed  to  allow  realistic  inclusion  of 
aerial  refueling  during  the  planning  phase,  so  that  MAC 
can  coordinate  with  the  Strategic  Air  Command  on  the  exact 
number  of  tankers  required  for  a  given  contingency.  This 
research  will  provide  such  a  tool. 

Present  Models  and  Concepts 

The  study  of  incorporating  aerial  refueling  into 
present  operations  plans  is  a  continuing  concern  to 
HQ  MAC/XP,  the  operations  planning  function  at  MAC  Head¬ 
quarters.  The  Command  still  lacks  definitive  guidance, 
however,  on  how  AR  should  be  optimally  used.  Three  cri¬ 
teria  have  bee-  i  used  in  the  past  when  considering  aerial 
refueling  for  mission  during  the  planning  phase: 

1.  First,  AR  is  considered  if  the  mission  must 
aerial  refuel  to  avoid  enroute  stops.  This  might  be  to 


decrease  overall  closure  time  or  the  fact  that  a  stop 
might  not  be  possible  due  to  fuel  availability  or  politi¬ 
cal  considerations.  In  any  case,  these  missions  receive 
first  priority  for  AR  during  the  planning  phase. 

2.  An  operational  or  aircraft  limitation  will 
place  a  mission  in  the  next  priority  for  AR.  For  example, 
either  runway  length  or  maximum  aircraft  gross  weight 
limitations  may  not  allow  the  airlifter  to  depart  with 
enough  fuel  to  reach  its  destination . 

3.  AR  is  sometimes  used  solely  to  decrease  the 
closure  time  of  an  operation.  This  criterion  is  not  used 
extensively  because  of  the  large  demand  such  an  operation 
may  have  on  tankers  from  other  users. 

These  practices  have  been  reasonably  effective  in  the  past; 
however,  an  optimization  method  of  selecting  missions  for 
AR  would  be  much  more  gseful  in  terms  of  minimizing  the 
use  of  tanker  resources  and  fuel  consumption. 

MAC  also  has  a  simulation  model  that  uses  the  AR 
mission  scenario  to  a  limited  extent.  This  model  is  called 
M-14  and  its  primary  purpose  is  to  test  the  capacity  of 
the  airlift  system  under  varying  conditions.  It  makes 
use  of  Monte  Carlo  simulation  techniques  to  model  all 
facets  of  MAC  airlift,  incorporating  all  strategic  airlift 
aircraft  and  including  such  factors  as  crews,  air  bases, 
materiel  depots,  and  logistics  materiel  used  in  daily  MAC 
operations  (Ref  17) .  M-14  does  include  the  AR  option  for 
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its  mission  scenarios;  however,  the  missions  chosen  for 
AR  are  done  so  randomly  with  no  attempt  at  optimizing  the 
use  of  tankers.  A  total  of  25  arbitrary  AR  rendezvous 
points  are  built  into  the  model,  and  the  model  selects  a 
rendezvous  point  for  a  specific  mission  not  to  minimize 
fuel  consumption,  but  to  minimize  the  fuel  upload  at  the 
cargo  onload  station.  Also,  tankers  are  not  modeled  in 
M-14,  but  instead  are  treated  as  unlimited  resources. 
Furthermore,  consideration  for  multiple  tanker  AR  is  not 
allowed  (Ref  17) .  A  model  which  selects  optimum  rendezvous 
points  and  tanker  departure  bases  is  needed  to  enhance 
M-14 1 s  treatment  of  the  AR  scenario  (Ref  13).  It  would 
provide  a  more  realistic  and  fuel  conservative  basis 
within  M-14  for  using  AR,  compared  to  the  present  arbitrary 
method  of  selecting  refueling  points  and  fuel  offloads. 

Research  results  reported  by  Captain  Bordelon  and 
Major  Marcotte  (Ref  5)  used  the  Simulation  Language  for 
Alternative  Modeling  (SLAM)  package  to  model  the  AR  mis¬ 
sion  scenario  and  were  successful  in  finding  an  optimum 
rendezvous  point  for  one  specific  scenario.  The  purpose 
of  their  research,  however,  was  to  define  a  general  AR 
policy  rather  than  develop  a  model  suitable  for  operational 
use.  Also,  optimization  techniques  were  not  used,  but 
rather  a  total  enumeration  of  feasible  rendezvous  points 
was  done  followed  by  the  selection  of  the  point  requiring 
the  least  amount  of  fuel.  Certain  assumptions  were  also 


made  that  limited  the  operational  usefulness  of  the  model. 
For  example,  multiple  tanker  AR  was  not  allowed.  Also,  a 
typical  scenario  was  built  into  the  model  with  no  provi¬ 
sion  for  input ing  actual  route  segments.  Although  unsuc¬ 
cessful  in  defining  a  general  AR  policy,  the  Bordelon- 
Marcotte  study  did  develop  a  valid  fuel  consumption  model, 
portions  of  which  are  used  in  this  study. 

Problem  Statement 

Airlift  planning  for  future  contingencies  must 
include  the  aerial  refueling  option.  Furthermore,  given 
limited  tanker  resources  and  fuel,  aerial  refueling  mission 
planning  must  be  done  as  efficiently  as  possible.  Liter¬ 
ally  every  new  aircraft  entering  the  Air  Force  inventory 
will  be  AR-capable;  therefore,  tanker  resources  will  con¬ 
tinue  to  be  limited.  Airlift  planners  presently  do  not 
have  a  method  available  to  them  that  will  provide  an 
optimal  solution  to  the  planned  rendezvous  point  and  tanker 
departure  base  problem.  A  model  that  can  do  this  will  not 
only  enable  airlift  planners  to  better  define  MAC  tanker 
requirements,  but  will  also  result  in  a  substantial 
overall  fuel  savings. 

Objectives  of  This  Study 

The  overall  objective  of  this  study  is  to  develop 
a  model  that  will,  given  a  route  of  flight,  enroute  wind 
factors,  airlifter  takeoff  fuel  load,  available  tanker 
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bases,  and  bases  usable  for  diversion  if  AR  is  unsuccess¬ 
ful,  determine  the  optimal  AR  rendezvous  point  and  tanker 
departure  base  to  minimize  the  total  fuel  consumed  by 
both  air lifter  and  tanker.  To  meet  the  objective,  the 
model  must  allow  the  user  to  select  any  combination  of  air¬ 
craft  (KC-10A,  C-5A,  C-141B,  KC-135A)  and  must  allow 
multiple  tanker  aerial  refueling.  Furthermore,  this  model 
must  be  capable  of  being  used  operationally  for  planning 
routine  AR  missions  as  well  as  for  planning  contingency 
airlift.  It  must  also  be  adaptable  for  use  with  MAC'S 
M-14  model  of  the  airlift  system. 

A  subobjective  of  this  study  involves  exercising 
the  model  in  accordance  with  the  experimental  design, 
resulting  in  an  analysis  of  the  sensitivity  of  total  fuel 
consumed  to  changes  in  airlifter  takeoff  fuel  and  cargo 
load.  This  analysis  serves  to  validate  the  results  of  the 
model  while  providing  the  planner  with  some  insight  as  to 
what  causal  factors  have  the  most  effect  on  minimizing 
total  fuel  consumed.  The  following  three  questions  are 
addressed  in  the  analysis. 

1.  Considering  airlifter  takeoff  fuel  and  air- 
lifter  cargo  load,  which  factor  has  the  most  effect  on 
total  fuel  consumed?  It  is  hypothesized  that  cargo  load 
will  have  the  most  effect. 

2.  Is  there  a  direct  causal  relationship  between 
takeoff  fuel  and  total  fuel  consumed?  MAC'S  present 
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policy  is  to  depart  the  airlifter  with  the  maximum  pos¬ 
sible  fuel  load  consistent  with  airfield  and  airlifter 
limitations.  This  practice  will  be  evaluated  as  to  its 
effects  on  minimizing  total  fuel  consumption  and  mini¬ 
mizing  tanker  support. 

3.  Using  a  normal  rendezvous  and  offload  method 
of  aerial  refueling  (as  opposed  to  the  "buddy"  rendezvous/ 
offload  operation  in  which  the  tanker  accompanies  the  air¬ 
lifter  over  a  portion  of  its  flight) ,  does  the  KC-10A 
offer  a  significant  fuel  savings  over  the  KC-135A?  It  is 
hypothesized  that  the  findings  will  show  the  KC-135A  more 
efficient  (in  terms  of  fuel  consumption)  during  normal 
rendezvous  AR. 

It  should  be  recognized  at  this  point  that  the 
objective  of  minimizing  total  fuel  consumed  may  not  always 
be  consistent  with  minimizing  the  use  of  tankers.  For 
this  reason,  tankers  are  modeled  as  resources  which  are 
available  as  needed,  but  which  add  a  penalty  (fuel  con¬ 
sumed)  for  each  additional  tanker  used.  The  model  balances 
the  effect  of  higher  airlifter  fuel  consumption  against 
additional  tanker  enroute  fuel  consumption  using  dynamic 
programming  techniques.  The  result  is  the  minimum  fuel 
consumed;  however,  the  number  of  tankers  used  in  support¬ 
ing  the  mission  may  not  be  the  minimum.  A  second  subobjec¬ 
tive  involves  addressing  this  question. 
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As  a  third  and  final  subobjective,  further  model 
validation  is  accomplished  by  comparing  model  derived 
optimal  rendezvous  points  and  tanker  departure  bases  with 
those  actually  used  on  operational  missions.  Model  output 
will  be  compared  with  actual  computer  flight  plan  data 
from  previous  AR  missions. 

Scope  and  Limitations 

This  study  is  limited  to  only  single-ship  air- 
lifter  cells;  however,  multiple  tanker  cell  refuelings 
will  be  allowed  at  each  possible  rendezvous  point.  The 
model  also  allows  only  one  type  of  tanker  (either  KC-10A 
or  KC-135A)  for  each  single-ship  airlifter  mission.  Also, 
"buddy"  AR  missions  are  not  included  in  this  model;  that 
is,  the  tankers  must  rendezvous  with  the  airlifter,  off¬ 
load  their  fuel,  then  return  to  their  respective  bases, 
insuring  minimum  overhead  destination  fuel  requirements 
are  met.  Both  C-5A  and  C-141B  aircraft  are  included  in 
the  model. 

The  model  requires  the  user  to  have  and  input  the 
following  data; 

1.  Airlifter  route  of  flight  using  geographic 
coordinates  (maximum  of  20  points) . 

2.  Tanker  departure  bases  to  be  considered,  also 
using  geographic  coordinates  (maximum  of  10  bases) . 
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3.  Airlifter  diversion  bases  to  be  used  in  case 
AR  is  unsuccessful  (maximum  of  10  bases  over  the  entire 
route  of  flight) . 

4.  Wind  factor  for  each  leg. 

5.  Type  of  airlifter  (C-5A  or  C-141B)  and  type  of 
tanker  (KC-135A  or  KC-10A) . 

6.  Airlifter  cargo  load,  takeoff  fuel  load,  and 
required  overhead  destination  fuel.  This  will  require  a 
small  amount  of  preliminary  fuel  planning  on  the  part  of 
the  user . 

The  model  then  determines  where  the  airlifter 
should  refuel,  the  tanker  offload,  and  the  tanker  departure 
base  to  be  used  in  order  to  minimize  total  fuel  consumed 
by  both  airlifter  and  tanker.  It  also  insures  that  the 
airlifter  has  enough  fuel  to  divert  to  the  closest  diver¬ 
sion  base  if  AR  is  unsuccessful. 

Benefits  Derived  From  This  Study 

The  justification  for  this  effort  falls  into  three 
general  areas.  First,  and  most  important,  it  will  help 
MAC  planners  to  define  a  Command  policy  on  aerial  refuel¬ 
ing,  leading  to  a  realistic  definition  of  tanker  require¬ 
ments.  Second,  it  can  later  be  incorporated  into  the  M-14 
computer  model  to  provide  a  sound  basis  for  which  roissjc.  ^ 
are  selected  for  AR.  And  third,  it  will  provide  an  inter¬ 
active  computer  tool  to  peacetime  operational  AR  planners 
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that  will  determine  optimal  rendezvous  points  and  tanker 
departure  bases.  This  will  result  in  more  training  per 
flight  hour,  and  should  result  in  considerable  fuel  savings. 
Also,  since  the  model  can  be  adapted  for  use  with  other 
aircraft,  the  results  of  this  study  can  be  used  in  pro¬ 
viding  a  basis  for  determining  future  tanker  force 
requirements . 

Overview 

Chapter  II  begins  with  a  brief  summary  of  normal 
airlift  fuel  planning  requirements,  followed  by  a  descrip¬ 
tion  of  the  additional  fuel  planning  required  by  aerial 
refueling.  The  aerial  refueling  scenario  is  then  concep¬ 
tualized  as  a  system  and  the  relationship  between  control 
and  response  variables  is  examined.  The  assumptions  made 
during  model  conceptualization  are  then  outlined  and 
justified. 

Chapter  III  deals  with  the  quantification  of  model 
parameters  and  the  computerization  of  the  model.  It  begins 
with  a  description  of  the  dynamic  programming  procedure 
as  applied  to  the  AR  scenario,  to  include  the  definition 
of  the  return  and  transformation  functions.  The  "Complex" 
Method  of  Box  optimization  algorithm  (Ref  6)  used  in  the 
model  is  then  described,  followed  by  the  derivation  of 
fuel  consumption  models  for  the  C-141B,  C-5A,  KC-135A,  and 
KC-10A.  A  discussion  of  validation  and  verification  com¬ 
pletes  the  model  description  in  Chapter  III. 


12 


Chapter  IV  presents  the  experimental  design  and 
results,  using  a  typical  airlift  scenario  for  the  analysis. 
No  attempt  is  made  to  derive  general  rules  to  optimize 
fuel  consumption,  but  rather  the  results  are  used  to 
verify  and  validate  the  model.  Model  output  is  explained 
and  compared  with  expected  results.  A  sensitivity  analysis 
of  certain  model  parameters  on  model  output  is  also  pre¬ 
sented.  Chapter  V  completes  the  study  with  conclusions 
and  suggestions  for  further  research. 


II  System  Description 

Before  describing  the  aerial  refueling  scenario 
as  a  system,  a  brief  overview  of  the  MAC  fuel  planning  pro¬ 
cess  is  required.  Detailed  fuel  planning  requirements  are 
outlined  in  the  respective  aircraft  Technical  Orders  and 
MAC  regulations  (Refs  7,  8,  18,  19).  This  discussion  will 
be  limited  to  the  details  required  for  justifying  the 
present  research  design.  Figure  1  graphically  portrays 
the  normal  fuel  planning  sequence  for  a  typical  airlift 
mission.  Figure  2  shows  the  additional  fuel  requirements 
for  an  AR  airlift  mission.  Figure  3  shows  a  portion  of 
an  actual  computer  flight  plan  depicting  the  different 
fuel  planning  phases.  These  phases,  numbered  from  one  to 
thirteen,  are  also  shotfn  on  Figure  1  for  reference. 

Normal  Fuel  Planning 

Phases  one,  two  and  three  of  the  normal  fuel  plan¬ 
ning  sequence  deals  with  the  enroute  portion  of  the  flight. 
Enroute  fuel  is  obtained  from  the  applicable  fuel  plan¬ 
ning  documents  (Ref  19)  based  on  aircraft  gross  weight, 
temperature  deviation  from  standard,  altitude  profile, 
true  airspeed,  and  flight  times.  Enroute  reserve  is 
added  to  account  for  unexpected  enroute  diversions  or 
inaccurate  wind  forecasts.  This  figure  is  always  defined 
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Normal  Fuel  Planning 


(Ref  18)  as  10  percent  of  the  fuel  required  over  the  por¬ 
tion  of  the  route  where  enroute  navigation  aids  are  inade¬ 
quate  to  accurately  position  the  aircraft  at  least  once 
each  hour  (also  known  as  the  Category  I  portion  of  the 
route) .  Enroute  fuel  also  includes  the  fuel  required  to 
climb  to  initial  level-off  altitude,  and  is  highly  sensi¬ 
tive  to  aircraft  gross  weight  (which  limits  altitude  capa¬ 
bility)  and  true  airspeed. 

Alternate  fuel,  phase  four,  includes  the  fuel 
required  to  climb  away  from  original  destination  and 
cruise  to  a  suitable  alternate.  The  altitude  and  true 
airspeed  used  in  planning  are  determined  by  the  distance 
from  original  destination  to  the  alternate.  This  fuel 
figure  is  again  obtained  from  the  applicable  fuel  planning 
documents  (Ref  19)  using  aircraft  gross  weight  overhead 
destination,  temperature  deviation  from  standard,  and  air 
distance  to  the  alternate  as  inputs.  An  alternate  is 
always  required  for  overseas  airlift  missions  (Ref  18) . 

Holding  fuel,  phase  five,  is  defined  as  the  fuel 
required  to  hold  for  45  minutes  at  10,000  feet  at  the 
alternate,  if  required,  or  to  hold  1  hour  and  15  minutes 
at  20,000  feet  at  the  destination.  This  fuel  figure  is 
based  on  aircraft  gross  weight  at  either  destination  or 
alternate  and  is  obtained  from  aircraft  fuel  planning  docu¬ 
ments  (Ref  19) . 
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Approach  and  landing  and  taxi/run-up  fuel,  phases 
six  and  nine,  are  assumed  constant  during  the  planning 
phase.  For  the  C-5A,  the  fuel  required  for  approach  and 
landing  is  5200  pounds,  and  the  fuel  required  for  taxi/ 
run-up  is  2800  pounds;  for  the  C-141B,  the  appropriate 
figures  are  2500  and  1900  pounds,  respectively. 

Phase  seven,  the  identified  extra  fuel,  includes 
the  fuel  added  for  any  other  reason  not  covered  in  the 
other  phases.  This  would  include  planned  additional  taxi 
time,  known  enroute  diversions  or  delays  due  to  traffic  or 
weather,  or  extra  fuel  carried  to  avoid  refueling  at  des¬ 
tination.  For  most  missions,  this  figure  is  zero. 

Phase  13,  the  required  overhead  destination  fuel, 
is  sometimes  referred  to  as  the  most  important  figure  in 
terms  of  enroute  fuel  management.  The  crew  must  periodi¬ 
cally  check  this  planned  figure  against  what  the  actual 
overhead  destination  fuel  will  be  throughout  the  flight. 
This  figure  includes  the  fuel  needed  to  fly  from  destina¬ 
tion  to  alternate,  hold,  and  approach  and  land  at  the 
alternate.  Therefore,  if  at  any  time  the  crew  determines 
that  the  actual  overhead  destination  fuel  will  be  less 
than  this  planned  figure,  it  must  take  some  action  to 
decrease  fuel  consumption  or  consider  diversion. 

Phases  8  and  10  are  totals,  and  are  shown  in 
Figure  1.  Phase  11  indicates  planned  actual  fuel  load. 
Phase  12  is  no  longer  included  on  computer  flight  plans. 


but  represents  the  difference  between  phases  10  and  11, 
or  unidentified  extra  fuel. 

AR  Mission  Scenario 

The  typical  AR  mission  scenario  and  fuel  require¬ 
ments  are  outlined  in  Figure  2.  Start,  taxi,  takeoff, 
enroute  reserve,  fuel  to  alternate,  holding,  and  approach 
and  landing  fuel  requirements  are  exactly  the  same  as 
those  for  a  normal  airlift  mission.  Cruise  fuel  planning 
for  an  AR  mission  normally  requires  several  entries  into 
the  enroute  fuel  charts  of  the  applicable  fuel  planning 
document  (Ref  19)  for  the  different  altitude  and  airspeed 
combinations  encountered  throughout  the  mission.  The  fol¬ 
lowing  discussion  will  illustrate  the  process  involved. 

Initial  climb  fuel  from  departure  point  is  normally 

computed'  to  optimum  altitude,  or  the  altitude  which  pro- 

« 

vides  a  300  foot-per-minute  rate  of  climb,  unless  AR  is 
anticipated  immediately  after  level-off.  If  this  is  the 
case,  level-off  is  accomplished  at  25,000-26,000  feet.  It 
has  been  shown  that  this  altitude  range  provides  the  best 
flight  control  response  for  both  airlifter  and  tanker  at 
typical  refueling  gross  weights.  The  refueling  track 
itself  is  defined  along  the  airlifter 's  route  of  flight, 
beginning  at  the  initial  point  (IP)  and  ending  at  the  end- 
AR  point.  The  control  point  (CP)  is  normally  100  nautical 
miles  from  the  IP  and  defines  the  geographical  point  the 
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Fig.  2.  AR  Mission  Fuel  Planning 
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Fig.  3.  Computer  Flight  Plan  Fuel  Block 


tanker  and  airlifter  ideally  should  have  completed  their 
rendezvous.  The  airlifter  crosses  the  IP  at  cruise  altitude 
and  airspeed  and  establishes  radio  contact  with  the  tanker, 
which  should  be  orbiting  at  the  CP.  At  80  nautical  miles 
from  the  CP,  the  airlifter  descends  to  1000  feet  below 
tanker  altitude  (unless  already  there) ,  and  the  tanker 
departs  the  orbit  and  begins  a  turn  towards  the  airlifter. 
The  aircraft  then  continue  flying  toward  each  other  until 
a  pre-established  separation  is  reached  (usually  8 
nautical  miles) ,  at  which  time  the  tanker  performs  a  180 
degree  turn  back  to  refueling  heading,  hopefully  rolling 
out  directly  in  front  of  the  airlifter  at  2  nautical  miles 
horizontal  and  1000  feet  vertical  separation.  Fuel  for 
this  entire  procedure  is  determined  by  entering  the  fuel 
planning  charts  with  a  planned  gross  weight,  AR  altitude, 
and  temperature  deviation  from  standard  at  altitude.  Air 
refueling  commences  at  the  CP  at  252  knots  indicated  air¬ 
speed  and  continues  to  the  end-AR  point.  The  end-AR  point 
is  typically  150-400  nautical  miles  from  the  CP,  depending 


on  the  number  of  tankers  in  the  refueling  cell.  After  AR 
is  complete,  the  airlifter  climbs  again  to  its  optimum 
altitude,  which  can  now  be  as  much  as  10,000  feet  lower 
than  its  pre-AR  optimum.  This  lower  altitude  capability, 
because  of  the  higher  airlifter  gross  weight,  becomes  a 
significant  factor  in  locating  the  optimum  refueling  points. 

One  additional  fuel  requirement  for  the  AR  mission 
scenario  requires  emphasis.  The  airlifter  must  always 
have  enough  fuel  at  the  end-AR  point  to  divert  to  a  suit¬ 
able  alternate  base  if  AR  is  unsuccessful  (Ref  18) .  This 
requirement  defines  a  lower  bound  on  the  fuel  state  of  the 
airlifter,  as  described  later  in  Chapter  III. 

Tanker  Fuel  Requirements 

Tanker  fuel  requirements  are  conceptually  the  same 
as  those  for  the  airlifter.  The  tanker  must  have  enough 
fuel  to  fly  to  the  CP,*orbit  for  one  hour  (in  case  the  air¬ 
lifter  arrives  late  at  the  IP) ,  return  to  its  recovery 
base,  hold  for  45  minutes,  approach,  and  land  with  a  spe¬ 
cific  fuel  reserve.  Table  1  outlines  these  requirements 
as  assumed  for  the  model.  These  figures  are  assumed  con¬ 
stant  for  each  tanker  mission,  leaving  the  enroute  cruise 
fuel  as  the  only  variable  fuel  figure.  The  enroute  cruise 
fuel  is  a  function  of  gross  weight  and  true  airspeed  and 
is  determined  in  the  same  manner  as  that  for  the  airlifter. 
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TABLE  1 


TANKER  RESERVE  REQUIREMENTS 


KC-10A 

KC-135A 

One  Hour  Orbit 

17322  lb 

11263  lb 

Climb,  Descent, 

Approach  and 

Landing 

13900  lb 

7800  lb 

45  Minutes 

Holding 

11500  lb 

6  500  lb 

Total 

42722  lb 

25563  lb 

System  Conceptualization 

The  aerial  refueling  scenario  is  modeled  in 
Figure  4  as  a  network,  with  each  node  representing  a  spe¬ 
cific  geographic  coordinate  along  the  airlifter's  route  of 
flight.  Node  1  then  becomes  overhead  destination,  and  node 
N  is  the  point  the  airlifter  assumes  its  cruise  altitude 
and  airspeed  profile  after  departure.  The  airlifter  enters 
the  network  with  a  specified  takeoff  fuel  and  must  arrive 
overhead  destination  with  at  least  planned  required  over¬ 
head  destination  fuel.  At  each  node,  the  airlifter  can 
either  refuel  or  choose  not  to  refuel;  if  the  AR  option  is 
chosen,  the  airlifter  must  then  decide  how  much  to  offload 
from  the  tanker.  Each  node  also  has  penalties  associated 
with  the  decisions  made,  in  terms  of  total  fuel  consumed. 
The  first  is  the  fuel  required  for  the  airlifter  to  fly 
the  next  segment  of  the  route  of  flight.  This  is  highly 
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dependent  on  the  refueling  decision,  because  of  the  strong 
direct  relationship  between  gross  weight  and  fuel  con¬ 
sumption.  The  second  penalty  is  the  fuel  consumed  by  the 
tanker  to  fly  from  its  departure  base  to  the  node.  This 
again  is  highly  dependent  on  the  refueling  decision.  The 
objective  of  the  model  is  to  determine  what  decisions 
should  be  made  at  each  node  to  minimize  the  total  of  these 
penalties . 

The  problem  is  further  modeled  as  a  system  in 
Figure  5,  showing  the  control  variables  and  their  effect 
on  the  response  variables.  The  system  boundaries  include 
the  airlifter  route  of  flight  and  all  available  tanker 
departure  bases  specified  as  input.  Tankers  are  available 
on  demand,  but  at  a  significant  cost  to  the  system.  The 
measures  of  effectiveness  of  the  system  are  the  response 
variables — total  fuel  consumed  and  number  of  tankers  used. 
The  control  variables,  for  a  given  scenario,  are  the  type 
of  tanker  used,  takeoff  fuel  weight,  cargo  load,  and  refuel¬ 
ing  decisions  possible  at  each  node.  The  variables  fixed 
for  a  given  scenario  are  the  type  of  airlifter,  airlifter 
route  of  flight,  wind  factor  for  each  leg,  airlifter  diver¬ 
sion  bases,  and  tanker  departure  bases. 

An  analysis  of  the  causal  loop  diagram  shows  the 
refueling  decision  to  be  the  driving  factor  for  all  other 
control  variables  affecting  the  total  fuel  consumed.  The 
decision  to  refuel  more  decreases  the  takeoff  fuel 
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requirement,  but  in  turn  increases  the  gross  weight  of  the 
airlifter,  thereby  indreasing  total  fuel  consumption. 

Also,  as  the  AR  decisions  increase,  the  number  of  tankers 
required  increases,  also  increasing  total  fuel  consumption. 
It  may  seem  that  cargo  load  has  an  independent  effect;  how¬ 
ever,  further  analysis  shows  that  increasing  cargo  weight 
increases  gross  weight,  which  in  turn  increases  the  refuel¬ 
ing  decisions  and  again  increasing  total  fuel  consumption. 
This  diagram  and  its  interrelationships  serve  to  point  out 
the  complexity  of  the  problem  to  be  solved. 

Simplifying  Assumptions 

Because  of  the  complex  interactions  between  the 
many  variables  in  the  problem,  certain  simplifying  assump¬ 
tions  were  made  and  must  be  justified  before  describing 
the  details  of  model  formulation.  These  assumptions  are 
as  follows: 

1.  Tankers  will  always  depart  with  maximum  fuel 
load.  This  assumption  is  consistent  with  alert  aircraft 
fuel  loads  and  the  objective  of  minimizing  tanker  use. 

The  tanker  can  offload  any  amount  of  fuel  to  the  airlifter, 
as  long  as  it  retains  enough  fuel  to  satisfy  reserve 
requirements  as  specified  earlier.  These  reserve  require¬ 
ments  are  assumed  constant  for  each  tanker  mission. 

2.  The  tanker  and  airlifter  fuel  for  climb, 
descent,  required  overhead  destination,  and  approach  and 
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landing  is  assumed  constant  during  the  optimization.  The 
only  variable  to  be  optimized,  then,  is  enroute  cruise 
fuel.  This  is  consistent  with  actual  mission  fuel  plan¬ 
ning  practices. 

3.  Tankers  will  fly  great  circle  direct  routings 
from  their  departure  bases  to  the  AR  control  point  and 
will  recover  at  the  same  base  in  the  same  manner.  The  air- 
lifter  also  will  fly  great  circle  routings  between  each 
point  on  its  route  of  flight. 

4.  The  possibility  of  AR  is  only  allowed  at  each 
geographical  point  defining  the  airlifter's  route  of 
flight.  Depending  on  the  length  of  the  mission,  this 
allows  a  decision  to  be  made  every  150-600  nautical  miles. 
Since  the  typical  aerial  refueling  track  is  250-500  miles 
long,  this  assumption  does  not  have  a  significant  effect 
on  the  optimal  solution.  It  does,  however,  considerably 
decrease  the  computational  time  required. 

5.  Fuel  used  by  the  airlifter  during  the  aerial 
refueling  procedure  is  assumed  to  be  the  same  as  if  the 
airlifter  were  at  enroute  cruise.  The  offsetting  effects 
of  lower  true  airspeed  during  the  refueling  maneuver 
(resulting  in  lower  fuel  consumption)  and  lower  altitude 
(resulting  in  higher  fuel  consumption)  make  this  assumption 
reasonable . 

6.  Wind  factors  are  assumed  to  remain  constant 
over  each  leg  of  the  airlifter's  route  of  flight.  This 
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assumption  is  again  consistent  with  normal  fuel  planning 
practices . 

These  assumptions,  together  with  the  limitations 
presented  in  Chapter  I  and  the  system  description  con¬ 
tained  in  this  chapter,  fully  define  the  problem  to  be 
solved.  In  summary,  the  problem  can  now  be  defined  as  a 
network,  with  nodes  representing  points  along  the  air- 
lifter's  route  of  flight,  and  the  output  of  the  system 
equal  to  or  greater  than  the  required  overhead  destina¬ 
tion  fuel.  Each  node  has  a  decision  associated  with  it 
and  a  cost  associated  with  that  decision.  The  objective 
is  to  minimize  the  total  cost  of  the  system. 


i 
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Ill  The  Model 


Dynamic  Programming  Concepts 

Dynamic  programming  is  a  mathematical  technique 
often  useful  for  making  a  sequence  of  interrelated  deci¬ 
sions.  It  provides  a  systematic  procedure  for  determining 
the  combination  of  decisions  that  maximizes  overall  effec¬ 
tiveness  (Ref  14:266).  Richard  Bellman  further  describes 
the  problem  ideally  suited  to  the  dynamic  programming  solu¬ 
tion  technique  (Ref  3:81): 

1.  In  each  case,  there  is  a  physical  system  charac¬ 
terized  at  every  stage  by  a  small  set  of  parameters,  the 
state  variables. 

2.  At  each  stage,  there  is  a  choice  of  a  number 
of  decisions. 

3.  The  effect  of  a  decision  is  a  transformation 
of  the  state  variables. 

4.  The  past  history  of  the  system  is  of  no  impor¬ 
tance  in  determining  future  actions. 

5.  The  purpose  of  the  process  is  to  maximize 
(minimize)  some  function  of  the  state  variables. 

Bellman's  fourth  point  is  the  most  significant  when 
applied  to  the  AR  optimization  problem.  It  says  that  what¬ 
ever  the  initial  decisions  are,  the  remaining  decisions 
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will  result  in  an  optimum  with  regard  to  the  first  decision. 
The  idea  of  sequencing  decisions  enables  a  complex  problem 
of  many  decisions  to  be  broken  down  into  a  sequence  of 
smaller  problems  with  only  a  few  variables  (Ref  16:183). 

Each  of  these  smaller  problems  is  called  a  state,  with  its 
optimal  decisions  resulting  in  a  return  to  the  system's 
total  objective  function.  Each  stage  also  has  a  state 
value  associated  with  it,  representing  the  state  of  the  sys¬ 
tem  at  that  stage.  The  stage  state  values  are  linked 
together  by  a  transition  function  which  defines  how  the 
system  state  changes  from  stage  to  stage.  Therefore,  the 
effect  of  a  policy  decision  at  each  stage  is  to  transform 
the  current  state  into  a  state  associated  with  the  next 
stage,  using  the  transition  function. 

Figure  6  shows  the  typical  dynamic  programming 
problem  in  a  serial  stage  format;  that  is,  each  stage 
occurs  in  direct  sequence.  Sn  represents  the  system  input 
state,  with  representing  the  state  value  at  stage  i. 

The  variable  d^  is  the  stage  decision,  and  r^  is  the  return 
at  stage  i  given  decision  d^.  is  the  state  transforma¬ 

tion  function,  and  SQ  the  system  output  state.  The  problem 
becomes  one  of  finding  the  optimum  of  a  function  of  the 
stage  returns  subject  to  constraints  on  the  decision 
variable: 
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Problem  Structure 


Minimize 


F(rl'r2 . rn> 

Subject  To  Gj  <  dj(S^)  £  H ^  j=l,2,...,n  (1) 


where  n  is  the  number  of  stages,  G^  is  the  lower  limit  of 
d j ,  the  decision  variable,  and  is  the  upper  limit.  The 
state  variables  are  in  turn  constrained,  either  between  con 
stants  or  functions  of  the  input  and  output  states. 


Formulation  of  the  AR  Problem 

The  AR  optimization  problem  is  ideally  suited  to 
the  use  of  dynamic  programming  solution  techniques.  To 
illustrate,  consider  a  problem  with  20  possible  AR  rendez¬ 
vous  points,  each  point  representing  a  decision  to  refuel 

or  not  to  refuel.  Using  network  analysis  optimization 

20 

techniques,  a  total  of  2  ,  or  1,048,576  possibilities 

would  have  to  be  considered  (Ref  2)  (see  Figure  7) .  Since 
dynamic  programming  effectively  eliminates  non-optimal  solu 
tions  at  each  stage  (recall  Bellman's  fourth  point  dis¬ 
cussed  in  the  previous  section) ,  only  a  small  subset  of 
the  total  number  of  possibilities  must  be  considered  at 
each  stage.  This  "principle  of  optimality"  allows  the  solu 
tion  of  a  much  larger  problem  with  considerably  less  com¬ 
puter  memory  and  time . 

Figure  6  can  again  be  used  to  illustrate  dynamic 
programming  formulation  of  the  AR  optimization  problem. 

Each  stage  represents  a  point  along  the  airlifter's  route 
of  flight,  and  the  state  variable  at  each  stage  is  the 
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airlifter' s  fuel  state  prior  to  making  the  refuel/no  refuel 
decision.  Therefore,  each  stage  is  also  a  possible  AR 
rendezvous  point.  The  decision  variable  represents  how 
many  tankers  are  needed  for  refueling  at  each  stage  based 
on  the  tanker's  available  fuel  for  offload.  For  example, 
a  decision  value  of  1.26  would  require  2  tankers — one 
tanker  would  offload  all  its  available  fuel  and  the  other 
only  .26  of  its  available  fuel.  A  decision  value  of  zero 
(the  model  actually  returns  a  value  of  1x10  ^ )  would  indi¬ 
cate  that  AR  would  not  take  place  at  that  point.  The  fuel 
available  for  offload  for  each  tanker  is  computed  by  sub¬ 
tracting  the  fuel  required  to  fly  from  departure  point  out 
to  the  AR  point  and  back  plus  required  tanker  reserves 
(see  Table  1),  from  the  tanker's  maximum  fuel  weight. 

The  return  function  for  each  stage  is  the  cost 
associated  with  the  decision  made  at  each  stage.  To  be 
consistent  with  the  objective  of  minimizing  total  fuel  con¬ 
sumed,  this  becomes  the  fuel  consumed  by  the  tanker  plus 
the  fuel  used  by  the  airlifter  to  fly  the  next  segment  to 
the  next  stage.  The  transition  function  transforms  the 
current  state  to  the  next  state  by  subtracting  the  fuel 
used  by  the  airlifter  between  the  stages  and  adding  the 
fuel  onload  from  the  tanker. 

At  each  stage,  the  airlifter  must  have  enough  fuel 
to  at  least  fly  the  next  route  segment  plus  divert  from  the 
next  stage  to  a  suitable  alternate  if  AR  is  unsuccessful 
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at  that  stage.  This  defines  the  lower  bound  on  the  deci¬ 
sion  value.  The  upper  constraint  is  driven  by  the  maximum 
allowable  gross  weight  or  maximum  fuel  load  of  the  air- 
lifter.  The  state  variables  are  constrained  by  the  same 
factors.  Further  details  of  how  the  constraints  are 
entered  in  the  model  are  located  in  Appendix  B. 

Quantification 

Aircraft  Parameters .  Table  2  lists  the  appropriate 
aircraft  limitations  and  parameters  used  in  the  model. 

These  figures  come  directly  from  the  respective  aircraft 
performance  manuals. 


TABLE  2 

AIRCRAFT  PARAMETERS 


Parameter 

C-5A 

C-141B 

KC-135 

KC-10A 

Max  Gross  Weight  (lb) 

712500 

323000 

270000 

590000 

Max  Fuel  Weight  (lb) 

318100 

153352 

163000 

34  5400 

Cruise  True  Airspeed  (kt) 

4  50 

432 

4  79 

479 

Operating  Weight  (lb) 

354000 

148800 

105000 

248470 

Fuel  Consumption  Models .  Gosch  and  Mooz  (Ref  12: 
20),  in  their  study  of  an  Air  Force  Fuel  Consumption  Model, 
develop  an  empirical  method  of  calculating  aircraft  fuel 
consumption  in  gallons  per  hour.  By  using  multiple  regres¬ 
sion  techniques,  they  were  able  to  derive  a  simple  relation 
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that  would  calculate  fuel  flow  as  a  function  of  aircraft 
gross  weight  and  true  airspeed.  The  general  form  of  the 
model  used  in  their  analysis  is 

FF  =  A  (GW)  B(TAS)  C  (2) 

where  A,  B,  and  C  are  constants,  GW  is  the  aircraft  gross 
weight  in  pounds,  TAS  is  the  true  airspeed  in  knots,  and 
FF  is  the  fuel  flow  in  gallons  per  hour.  Aircraft  alti¬ 
tude,  a  significant  factor  in  determining  fuel  consumption, 
is  accounted  for  by  the  gross  weight  term,  since  a  gross 
weight  defines  a  specific  altitude  capability.  This 
general  model  was  used  to  derive  specific  models  for  the 
C-5A  and  C-141B,  KC-135A,  and  KC-10A,  using  the  Statistical 
Package  for  the  Social  Sciences  (SPSS)  regression  routines 
(Ref  4).  The  resulting  models  for  each  aircraft  are: 


KC-10A: 

FF  =  7 . 995xl0~6 (GW) ' 9795 (TAS) 1 ' 1214 

(3) 

KC-135A: 

FF  = 2 . 54xl0~8 (GW) '916 (TAS) 2' 24 

(4) 

C-5A  and  C-141B: 

FF  =  1 . 81xl0-3 (GW) *  834 5 (TAS)  '5503 

(5) 

A  factor  of  6.5  pounds  per  gallon  is  used  to  convert  the 
fuel  flow  in  each  case  to  the  more  commonly  used  pounds  per 
hour.  A  full  description  of  the  regression  analysis  is 
presented  in  Appendix  A. 

The  fuel  actually  consumed  over  a  route  segment 
is  calculated  using  the  following  simple  relation: 


Fuel  Consumed  =  [XRT/ (TAS+WF) ]  *  FF,  (6) 

where  FF  is  the  fuel  flow  as  defined  above,  TAS  is  the  true 
airspeed  in  knots,  XRT  is  the  length  of  the  route  segment 
in  nautical  miles,  and  WF  is  the  wind  factor  over  that 
route  segment.  Appendix  B  describes  the  model  used  to 
calculate  the  great  circle  distance  between  two  points 
from  the  geographic  coordinates. 

Return  Function .  Using  the  analysis  in  the  pre¬ 
vious  section,  the  system  return  for  each  stage  is  the  fuel 
required  by  the  tanker  plus  the  fuel  required  by  the  air- 
lifter  to  fly  the  route  segment  from  that  stage,  or 

r±  =  <d  .  >  *  (XMAXFT-FTNKi)  +  [  (FF^*XRT^)/  (TAS+WF^)  ]  *6.5 

3  1  (7) 

where 

d^  =  decision  at  stage  i 

<d.>  =  integer  number  of  tankers  required  to  deliver 
1  the  offload.  <d.>isthe  integer  greater  than 
or  equal  to  dj  3 

XMAXFT  =  maximum  fuel  weight  (takeoff  fuel)  of  the 
tanker 

FTNK^  =  the  fuel  available  for  offload  at  stage  i 
r^  =  the  return  for  stage  i 


The  fuel  flow  is  calculated  using  the  airlifter's  gross 
weight  at  each  stage  after  refueling. 


Transition  Function.  The  state  transformation 


between  stages  is  determined  from  the  following: 

SMAXi  =  SMAXi+1  -  [ {FFi*XRTi) / (TAS+WFi) ] 

*  6.5  +  £di+1*FTNKi+1 J  (8) 

where  SMAX^  is  the  airlifter  fuel  state  at  stage  i  and 
SMAXi+1  is  the  airlifter  fuel  state  at  the  previous  stage. 
The  fuel  flow  is  again  calculated  using  the  airlifter' s 
gross  weight  at  each  stage  after  refueling. 

System  Return  Function .  The  system  return  func¬ 
tion  used  in  the  optimization  procedure  is  simply  the  sum 
of  the  individual  stage  returns  as  defined  in  eq  (7) ,  or 

n 

F  =  l  r.  (9) 

i=l  1 

• 

for  the  n  stage  problem. 

Tanker  Fuel  Available  for  Offload.  Fuel  available 
for  offload  to  the  airlifter  is  calculated  at  each  stage 
from  the  following  expression: 

FTNKj,  =  XMAXFT  -  { [  (FF*XTNK^)  /TAST]  *  2  +  RES}  (10) 

where 

XMAXFT  =  the  maximum  allowable  fuel  load  for  the 
tanker 

FF  =  tanker  fuel  flow  in  pounds  per  hour 


XTNK .  =  great  circle  distance  in  nautical  miles 
1  from  the  closest  tanker  departure  base  to 
stage  i 

TAST  =  tanker  true  airspeed  in  knots 
RES  =  required  tanker  reserves 

The  gross  weight  used  in  the  fuel  flow  calculation  is  the 
tanker  operating  weight  plus  one-half  of  the  maximum  fuel 
load.  This  results  in  an  average  fuel  flow  and  assumes 
that  the  tanker  burns  or  offloads  all  fuel  during  its  mis¬ 
sion. 

Solution  Algorithm 

The  dynamic  programming  algorithm  proceeds  as  fol¬ 
lows  (refer  to  Figure  6) : 

1.  At  stage  1,  the  upper  and  lower  limits  on  the 
input  fuel  state  are  evaluated.  The  interval  from  the 
high  limit  to  the  low  limit  is  divided  into  a  finite  series 
of  equidistant  steps. 

2.  Beginning  at  the  lower  limit  on  the  state,  S^ , 
the  optimum  return  is  found  over  the  range  of  decisions  d^ : 

opt 

fl(Sl)  =  dl  <r1(S1,d1) }  (11) 

The  "Complex"  Method  of  Box  (see  next  section  and  Appen¬ 
dix  D)  is  used  to  find  the  optimum  return  for  each  stage. 

3.  The  state  value  is  then  incremented  and  the 
optimum  is  found  for  this  value  of  S^. 
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4.  The  process  is  repeated  until  the  upper  bound 


on  is  reached. 

5.  Starting  with  n=2,  the  process  moves  backwards 
through  the  stage  sequence.  The  limits  on  Sn  are  evaluated 
and  the  state  interval  is  divided  into  equidistant  steps. 

6.  The  optimum  return  is  found  for  each  value  of 
S  by  the  formula 

opt 

f„,S)  “  d„{r„<W  +  £n-l(Sn-dn)}-  <12> 

7.  Recursion  is  followed  until  stage  n  is  reached 

and  the  optimum  found  for  all  values  of  S  . 

n 

8.  The  algorithm  then  follows  the  optimum  path 

forward  to  recover  the  optimum  decisions  as  follows:  The 

input  fuel  state  S  is  located  in  the  -able  of  state  values 

n 

for  stage  n,  with  the  corresponding  f  and  optimum  decision. 

♦ 

The  transition  function  between  S  and  S  is  used  to  find 

n  n-1 

the  optimum  Sn_^  and  its  corresponding  values.  This 
forward  path  is  followed  through  all  stages  until  and 
the  corresponding  r^  and  optimum  decisions  are  found  (Ref 
16:185)  . 

The  "Complex"  Method  of  Box 

The  "Complex"  Method  of  M.  J.  Box  (Ref  6)  finds  the 
optimum  value  of  a  multivariable,  nonlinear  function  sub¬ 
ject  to  nonlinear  inequality  constraints: 
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Min imize 


F  (X1#X2, . . . ,Xn) 

Subject  To  Gk  <  Xk  <  Hk  k=l,2,...,m  (13) 

The  implicit  variables  Xn+i . Xm  are  dependent  functions 

of  the  explicit  independent  variables  X^ ,X^ , . .  .  ,Xn .  The 
upper  and  lower  limits  Hk  and  Gk  are  either  constants  or 
functions  of  the  independent  variables.  The  method  itself 
is  a  sequential  search  technique  which  has  been  proven 
successful  in  solving  nonlinear  objective  functions  with 
nonlinear  inequality  constraints.  The  procedure  should 
tend  to  find  the  global  maximum  by  selecting  an  initial  set 
of  points  that  are  randomly  scattered  through  the  feasible 
region.  The  algorithm  evaluates  the  objective  function 
and  searches  through  the  feasible  region  for  an  optimum 
value,  starting  with  a  random  initial  set  of  points,  until 
it  satisfies  user  defined  convergence  criteria.  A  detailed 
description,  along  with  selected  solution  parameters,  is 
presented  in  Appendix  D. 

Computerization 

Appendix  C  contains  a  complete  FORTRAN  V  computer 
listing  of  the  AR  optimization  model,  to  include  explana¬ 
tions  of  all  variables  and  subroutines.  Sample  raw  and 
formated  output  is  presented  in  Chapter  IV  and  Appendix  E. 
User  input  instructions  are  included  as  comments  in  the 
main  DYNAM  program. 


41 


Verification  and  Validation 

Fishman  and  Kiviat  (Ref  11)  divide  the  process  of 
model  evaluation  into  three  categories: 

1.  Verification,  to  insure  that  the  model  behaves 
as  intended; 

2.  Validation,  to  test  agreement  between  the 
behavior  of  the  model  and  the  system  it  is  meant  to  por¬ 
tray  ;  and 

3.  Problem  analysis,  which  analyzes  and  interprets 
the  data  generated  by  the  experimental  design. 

In  other  words,  the  primary  concern  is  with  the 
internal  consistency  of  the  model,  its  correspondence  with 
the  real  system,  and  correct  interpretation  of  the  result¬ 
ing  data  (Ref  21:210).  History  has  seen  three  different 
viewpoints  emerge  on  the  correct  method  to  conduct  scien¬ 
tific  inquiry — the  rationalist,  the  empiricist,  and  the 
absolute  pragmatist  points  of  view  (Ref  21:212). 

The  rationalist  can  be  seen  on  one  end  of  the 
philosophical  spectrum  of  model  evaluation.  According  to 
the  rationalist,  the  validity  of  the  model  is  dependent 
solely  on  whether  one  accepts  the  basic  premises  and  logic 
used  to  construct  the  model.  No  proof  or  experimental  mea¬ 
surement  of  the  premises  is  required.  Because  of  this  lack 
of  proof,  the  rationalist  viewpoint  is  unacceptable  to  the 
empiricist  at  the  other  end  of  the  spectrum.  Empiricism 
refuses  to  allow  any  premises  or  assumptions  that  are  not 
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verified  by  experiment  or  analysis  of  experimental  data. 
This  way  of  thinking  is  extremely  limited  in  application, 
since  it  requires  a  model  to  begin  with  proven  or  verifi¬ 
able  facts,  and  assumptions  are  not  allowed  (Ref  21:210). 
The  third  theory  of  model  evaluation,  absolute  pragmatism, 
looks  at  the  model  as  a  black  box  and  is  only  concerned 
with  the  usefulness  of  the  output  of  the  black  box,  not  its 
internal  structure.  The  model,  then,  is  valid  only  if  it 
fulfills  the  purpose  for  which  it  was  designed.  Again, 
this  viewpoint  overlooks  the  very  important  conceptualiza¬ 
tion  phase  of  model  formulation.  What  is  needed  is  an 
approach  that  combines  the  merits  of  all  three  of  these 
philosophies,  a  sort  of  verification/validation  systems 
approach,  allowing  these  viewpoints  to  merge  and  interact. 
Shannon  (Ref  21:215)  calls  this  the  utilitarian  approach 
to  the  model  verification  process.  This  approach  is  used 
to  verify  and  validate  the  AR  optimization  model. 

The  first  stage  of  the  utilitarian  approach  is  to 
seek  face  validity  of  the  internal  structure  of  the  model 
based  on  a  priori  knowledge,  past  research,  and  existing 
theory  (Ref  21:215).  The  internal  structure,  interaction 
between  variables,  and  analysis  of  the  AR  scenario  as  a 
system  have  been  examined  in  the  conceptualization  phase 
with  a  view  towards  showing  face  validity  of  the  model. 
Assumptions  have  been  stated  and  justified.  The 
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application  of  dynamic  programming  and  the  "Complex"  Method 
of  Box  will  complete  this  stage  in  Chapter  IV. 

The  second  stage  is  also  concerned  with  the  valida¬ 
tion  of  the  internal  structure  of  the  model,  and  consists 
of  empirically  testing  the  assumptions  and  hypotheses  used 
(Ref  21:216).  This  stage  will  also  be  satisfied  in  the 
next  chapter  when  the  model  is  exercised  in  accordance  with 
the  experimental  design  and  model  output  is  compared  with 
expected  results.  Also,  certain  model  solution  parameters 
will  be  varied  to  test  the  model's  sensitivity  to  these 
parameters . 

The  third  stage  attempts  to  verify  the  model's 
ability  to  predict  the  behavior  of  the  real  world  system 
(Ref  21:216).  This  will  be  accomplished  by  comparing  model 
derived  optimal  rendezvous  points  and  fuel  consumption 
with  those  actually  used  on  operational  missions.  The 
philosophical  question  of  how  important  it  is  that  the 
model  be  a  true  and  isomorphic  reflection  of  the  real  sys¬ 
tem  will  also  be  addressed. 
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IV  Experimental  Design  and  Analysis 


The  purpose  of  this  chapter  is  to  outline  the 
design  of  the  experiment  used  to  validate  the  AR  optimiza¬ 
tion  model  and  to  present  the  results  of  the  analysis. 

This  will  be  done  in  four  parts.  First,  the  experiment 
used  to  analyze  the  effects  of  the  control  variables  on 
total  fuel  consumed  will  be  described,  followed  by  a  dis¬ 
cussion  of  the  three  questions  posed  in  the  first  subobjec' 
tive.  Then,  the  question  of  whether  minimizing  total  fuel 
consumed  is  consistent  with  minimizing  the  use  of  tankers 
will  be  addressed.  A  sensitivity  analysis  of  optimiza¬ 
tion  parameters  will  then  be  presented,  followed  by  a 
comparison  of  model  output  with  actual  operational  mission 
computer  flight  plans.*  The  combination  of  these  four 
analyses  will  complete  the  model  verification  and  valida¬ 
tion  as  outlined  in  Chapter  III. 

Part  One 

Experimental  Design .  A  typical  aerial  refueling 
scenario  was  chosen  for  analysis,  although  the  model 
accepts  any  possible  combination  of  airlifter  routings 
and  tanker  departure  bases.  Although  the  selection  of 
one  specific  routing  and  set  of  tanker  bases  limits  the 
results  of  this  analysis  in  terms  of  deriving  a  general 
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optimal  rule  for  selecting  AR  rendezvous  points,  it  does 
serve  to  validate  the  model  and  its  assumptions.  It  is 
felt  that  because  of  the  limitless  number  of  combinations 
of  routings  and  possible  tanker  departure  bases,  a  meaning¬ 
ful  general  rule  would  not  apply  in  all  cases. 

The  scenario  chosen  is  shown  in  Table  3  and 

involves  an  airlift  mission  from  Biggs  Army  Airfield, 

El  Paso,  Texas,  nonstop  to  Dhahran,  Saudi  Arabia.  The 

routing  was  taken  from  an  actual  computer  flight  plan  and 

reflects  International  Civil  Aviation  Organization  (ICAO) 

preferred  routings  through  Europe  and  the  Middle  East. 

Table  3  also  shows  the  assumed  wind  factors  (also  obtained 

from  the  computer  flight  plan)  and  the  model  generated 

great  circle  distances  between  each  point.  The  coordinates 

are  given  in  degrees-tenths  of  degrees  for  the  latitude 

and  longitude,  where  a  negative  coordinate  represents  east 
'  * 

longitude  and  south  latitude.  Table  4  shows  the  tanker 
departure  bases  and  airlifter  diversion  bases  used  in  the 
experiment.  Figure  8  shows  the  route  of  flight,  tanker 
departure  bases,  and  diversion  bases  geographically. 

Table  5  summarizes  the  experimental  design  factors 
and  levels  for  the  first  phase  of  analysis.  A  maximum 
takeoff  gross  weight  of  650,000  pounds  is  assumed  for  the 
C-5A  because  of  airfield  limitations;  for  the  C-141B,  the 
maximum  takeoff  gross  weight  is  assumed  to  be  318,800 
pounds.  Cargo  weights  for  both  aircraft  are  typical 
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TABLE  4 

EXPERIMENTAL  DESIGN  DATA 
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Airlifter  Route  of  Flight 


TABLE  5 


EXPERIMENTAL  DESIGN  LEVELS 


Factor 


Level 


Aircraft  Combination 


Airlifter  Cargo  Load 


Airlifter  Takeoff  Fuel 


a)  C-5A/KC-135A 

b)  C-5A/KC-10A 

C)  C-141B/KC-135A 

C-5A: 


a) 

100 

,000 

lb 

b) 

125 

,000 

lb 

c) 

150 

,000 

lb 

C- 

141B: 

a) 

20 

,000 

lb 

b) 

35 

,000 

lb 

c) 

50 

,000 

lb 

C- 

5A: 

a) 

46 

,000 

lb 

b) 

71 

,000 

lb 

c) 

96 
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lb 

d) 
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lb 

e) 
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f) 
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lb 

g) 
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lb 

g) 
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lb 

productive  loads  for  a  mission  of  this  type.  Takeoff  fuel 
figures  start  at  46,000  pounds  for  the  C-5A  and  30,000 
pounds  for  the  C-141B,  and  are  incremented  25,000  pounds 
and  15,000  pounds,  respectfully,  up  to  a  maximum  consistent 
with  the  selected  takeoff  gross  weights.  The  C-5A  is 
allowed  to  refuel  from  either  type  tanker;  however,  the 
C-141B  only  refuels  from  the  KC-135A.  The  C-141B/KC-10A 
combination  is  considered  infeasible  for  the  rendezvous/ 
offload  method  of  AR  because  of  the  tanker's  large  refuel¬ 
ing  capability  and  fuel  consumption,  and  the  airlifter's 
limited  ability  to  accept  large  amounts  of  fuel  at  one 
time.  Since  the  model  is  deterministic,  only  one  run  of 
each  of  the  63  data  points  was  made. 

Analysis.  Appendix  G  contains  a  listing  of  all 
model  output  for  the  experimental  design.  Figures  9,  10, 
and  11  contain  plots  of  total  fuel  consumed  and  fuel  off¬ 
loaded  for  the  C-5A/KC-135A  and  C-5A/KC-10A  combinations, 
all  as  a  function  of  takeoff  fuel  and  cargo  load.  Figures 
12,  13,  and  14  contain  the  same  plots  for  the  C-141B/KC-135A 
combination. 

For  all  air lif ter/tanker  combinations,  Figures  9 
through  14  show  that  cargo  load  appears  to  have  the  great¬ 
est  effect  on  total  fuel  consumed.  This  conclusion  was 
tested  statistically  for  the  C-5A  data,  the  results  of 
which  are  contained  in  Appendix  G.  As  was  expected,  cargo 
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Fig.  11.  C-5  Fuel  Consumed/Offloaded 

Cargo  Load  =  150,000  lbs 


Fig.  12.  C-141  Fuel  Consumed/Offloaded 

Cargo  Load  =  20,000  lbs 


Fig.  13.  C-141  Fuel  Consumed/Off loaded 

Cargo  Load  =  35,000  lbs 


load  had  a  much  greater  effect  on  the  total  fuel  consumed 
than  takeoff  fuel.  This  result  is  due  to  the  sensitivity 
of  the  airlifter's  fuel  consumption  on  total  gross  weight 
and,  hence,  cargo  load.  Changing  the  takeoff  fuel  weight 
many  times  just  changes  the  amount  of  fuel  offloaded  at 
one  or  more  aerial  refueling  points,  which  results  in  very 
small  changes  in  total  fuel  consumed.  Furthermore,  the 
points  chosen  for  AR  also  were  not  affected  by  takeoff 
fuel;  the  primary  factor  affecting  the  AR  decision  appears 
to  be  the  distance  of  the  nearest  tanker  base  from  the 
stage,  which  in  turn  affects  the  tanker's  fuel  consumption. 
Therefore,  a  stage  with  a  tanker  base  nearby  would  have  a 
much  better  chance  of  being  used  as  an  AR  point  than  one 
with  a  tanker  base  far  away. 

A  direct  causal  relationship  between  takeoff  fuel 
and  total  fuel  consumed  does  not  appear  to  exist.  In  most 
Of  the  combinations,  total  fuel  consumed  decreases  as  take¬ 
off  fuel  increases,  as  might  be  expected.  This  is  a  result 
of  the  fact  that  more  takeoff  fuel  in  most  cases  decreases 
the  need  for  AR.  Figures  9,  11,  12,  and  13  show  cases 
where  the  total  fuel  consumed  remains  constant  over  a  por¬ 
tion  of  the  takeoff  fuel  range.  This  result  is  due  to  the 
model  varying  the  fuel  offload  at  the  same  AR  point  without 
using  additional  tankers.  For  example,  in  Figure  9,  a 
takeoff  fuel  of  46,000  pounds  results  in  a  KC-10A  fuel 
offload  of  212,500  poinds  at  stage  1;  as  the  takeoff  fuel 
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increases  to  196,000  pounds,  the  KC-10A  offload  decreases 
at  stage  1  to  a  figure  of  62,500  pounds.  Since  the  aircraft 
gross  weight  remains  the  same  in  all  cases,  the  airlifter' s 
fuel  consumption  also  does  not  change.  The  tanker's  fuel 
consumption  also  stays  the  same,  since  the  refueling 
requires  one  tanker  in  each  case.  Therefore,  total  fuel 
consumption  remains  constant.  Figures  10,  12,  and  14  show 
instances  where  the  total  fuel  consumed  increases  as  take¬ 
off  fuel  increases.  This  can  be  explained  again  by  the  fact 
that  higher  gross  weights  require  more  fuel;  in  these  cases, 
the  additional  gross  weight  caused  by  the  higher  takeoff 
fuel  increased  fuel  consumption  to  the  point  that  an  addi¬ 
tional  aerial  refueling  was  required.  This  increased  the 
response  variable,  total  fuel  consumed. 

For  the  C-5A,  Figures  9  through  11  graphically  show 
the  KC-135A  to  be  more  efficient  than  the  KC-10A  for  the 
normal  rendezvous  mode  of  aerial  refueling.  This  is 
because  of  the  significantly  higher  fuel  consumption  of  the 
KC-10A  and  the  overall  inefficiency  of  offloading  all  of 
the  KC-lOA's  available  fuel  once  it  arrives  at  the  AR  point. 
Recent  operational  tests  show  the  KC-10A  most  efficient  in 
the  "buddy"  AR  method,  where  the  tanker  accompanies  the  air- 
lifter  over  its  route  of  flight,  refueling  the  airlifter 
as  required.  Using  the  KC-10A  as  a  refueler/airlifter  in 
this  manner  results  in  significant  fuel  savings  by 
reducing  airlifter  requirements  (by  carrying  up  to  three 


C-141B  or  one  C-5A  cargo  equivalent  loads)  and  simultane¬ 
ously  reducing  tanker  requirements  (by  accomplishing  all 
AR  requirements  along  the  way)  (Ref  9) . 

Figures  9  through  14  also  show  the  total  fuel  off¬ 
loaded  for  various  cargo  loads  and  takeoff  fuel  weights. 

As  is  expected,  fuel  offload  requirements  decreased  in  most 
cases  as  takeoff  fuel  increased.  Figure  12  shows  one 
exception,  where  the  offload  requirement  jumped  to  a  higher 
value  as  takeoff  fuel  increased  from  75,000  pounds  to 
90,000  pounds.  This  is  one  case  where  the  higher  gross 
weight  (hence  higher  fuel  consumption)  caused  an  additional 
AR.  In  all  other  cases,  the  higher  takeoff  fuel  effectively 
decreased  AR  requirements  along  the  route  of  flight.  For 
the  C-5A,  fuel  offload  requirements  were  always  a  little 
higher  for  the  KC-10A  than  for  the  KC-135A.  This  is 
indirectly  due  to  the  higher  fuel  consumption  of  the  KC-10A. 
Since  more  fuel  was  required  to  get  the  KC-10A  to  the  AR 
point,  this  additional  cost  to  the  system  forced  higher  off¬ 
loads  whenever  the  AR  option  was  selected.  These  higher 
offloads  increased  airlifter  fuel  consumption,  which  in  turn 
increased  offload  requirements  later  on. 

Part  Two 

The  objective  of  minimizing  fuel  consumption  is 
clearly  not  always  compatible  with  that  of  minimizing  the 
use  of  tankers.  Appendix  G  contains  the  number  of  tankers 


used  for  each  of  the  aircraft  combinations  in  the  experi¬ 
mental  design.  The  use  of  tankers  was  minimized  only  to 
the  extent  of  minimizing  total  fuel  consumed.  In  many 
cases,  the  mission  could  be  accomplished  with  fewer  tankers 
by  allowing  the  airlifter  to  offload  more  fuel  from  each 
tanker,  up  to  a  maximum  of  FTNK.  This,  of  course,  increases 
the  fuel  consumed  by  the  airlifter  and,  hence,  by  the  total 
system. 

From  the  results  of  the  data  used  in  the  analysis, 
the  real  benefits  derived  from  the  use  of  the  AR  optimiza¬ 
tion  model  are  modified  slightly  into  two  main  areas. 

First,  as  mentioned  in  Chapter  I,  its  use  during  peacetime 
operational  mission  planning  will  result  in  a  considerable 
fuel  savings  where  the  number  of  tankers  used  is  not  a 
critical  factor.  In  fact,  using  more  tankers  increases 
peacetime  training  missions  without  additional  fuel  costs. 
This  is  important  to  the  Strategic  Air  Command,  where  air¬ 
crew  members  have  traditionally  been  shorted  flying  hours 
due  to  fuel  costs.  Second,  the  model  can  be  used  during 
planning  for  large  scale  airlift  operations,  where  a  number 
of  airlifters  are  enroute  to  a  common  destination  via  the 
same  route  (these  are  known  as  airlifter  enroute  cells) . 
Unused  tanker  offload  capability  after  one  aerial  refuel¬ 
ing  can  be  used  to  refuel  a  subsequent  airlifter.  This 
type  of  operation  would  bring  together  the  two  objectives 
of  minimizing  tanker  use  and  fuel  consumed.  This  type  of 


■l..:,  .1 A. 


61 


analysis  is  also  what  is  needed  to  determine  tanker  force 
requirements  for  a  given  contingency. 


Part  Three 

Appendix  G  presents  the  scenario  used  and  the  data 
obtained  by  varying  the  number  of  iterations  used  in  the 
optimization  (ITMAX)  and  the  fuel  state  step  size  (NSTEP) . 
It  was  thought  that  these  two  parameters  had  the  most 
effect  on  the  efficiency  of  the  Box'  algorithm.  Other 
model  parameters  were  left  at  the  values  recommended  by 
Box  (Ref  6) . 

Varying  ITMAX  had  no  effect  on  the  results  of  the 
model  output,  at  least  for  this  scenario.  As  might  be 
expected,  the  convergence  criteria  stopped  the  iteration 
process  as  soon  as  the  complex  points  functional  values 
remained  within  .01  of  each  other  for  five  consecutive 
iterations.  When  ITMAX  was  increased  beyond  100,  the  com¬ 
puter  execution  time  remained  constant,  as  did  the  total 
fuel  consumed  figure.  This  indicates  that  the  algorithm 
was  behaving  as  intended. 

As  NSTEP  was  decreased,  the  step  size  between 
allowable  fuel  states  increased,  and  the  optimal  solution 
for  total  fuel  consumed  changed  slightly.  Decision  values 
only  changed  in  the  third  decimal  place,  and  the  optimal 
aerial  refueling  points  stayed  the  same.  The  changes  in 
fuel  offload  and  total  fuel  consumed  are  small  as  compared 
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to  the  actual  figure.  As  might  be  expected,  the  parameter 
NSTEP  only  tends  to  refine  the  model  output  and  can  be 
varied  by  the  user  as  needed  to  decrease  computer  execution 
time . 

Small  values  of  NSTEP  (less  than  50)  can  generate 
errors  of  2  percent  or  greater  between  consecutive  stage 
fuel  states.  This  is  because,  when  calling  back  the  opti¬ 
mal  solution  after  recursion,  the  model  tries  to  match 
state  values  between  stages  to  recover  the  next  stage 
optimal  solution.  If  it  cannot  match  the  state  values 
exactly,  it  selects  the  next  higher  state  value  and  its 
optimal  solution,  which  is  at  most  1/NSTEP  in  error. 
Therefore,  NSTEP  = 100  reduces  the  error  to  1  percent,  and 
this  value  was  selected  for  all  data  runs  in  this  study. 

In  general,  the  higher  NSTEP  is,  the  more  accurate  the 
model  output  will  be;  however,  computer  execution  time 
increases  in  direct  proportion  to  NSTEP. 

On  two  data  runs  with  NSTEP  equal  to  100,  the  model 
forced  an  AR  with  a  small  offload  of  less  than  the  state 
step  size.  This  was  also  caused  by  the  model's  inability 
to  exactly  match  state  values  during  the  optimal  solution 
call-back.  This  situation  can  be  corrected  by  increasing 
NSTEP  to  150 — the  model  then  combines  this  small  offload 
with  another  AR  later  along  the  route  of  flight.  In  both 
cases,  the  optimal  solution  was  not  affected  by  more  than 
1  percent. 
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Part  Four 

Six  typical  airlift  missions  were  selected  for  use 
in  comparing  model  output  with  actual  computer  flight  plan 
data.  These  missions,  input  data,  model  output,  computer 
flight  plan  information,  and  comments  on  each  comparison 
are  presented  in  Appendix  G.  This  section  contains  general 
*  comments  on  the  comparisons. 

Present  computer  flight  planning  procedures  appear 
extremely  effective  in  using  optimal  refueling  points, 
since  model  selected  AR  points  generally  agree  with  those 
used  on  actual  missions.  From  the  six  examples  selected 
for  analysis,  it  appears  that  fuel  is  wasted  by  offloading 
too  much  fuel  at  each  AR  point.  This  results  in  higher 
fuel  consumption  (approximately  25  percent  of  extra  fuel 
carried  is  burned  just  by  carrying  it)  and  excessive  fuel 
loads  overhead  destination.  Model  fuel  states  are  in 
close  agreement  with  those  on  the  computer  flight  plans, 
verifying  the  fuel  consumption  equations  used  in  the  model. 
The  small  differences  here  can  be  attributed  to  one  or  both 
of  the  following: 

1.  The  model  assumes  great  circle  distances 
between  stage  points.  The  computer  flight  plan  uses  airway 
or  route  distances  between  navigation  aids  as  published  on 
ICAO  aeronautical  charts.  The  computer  flight  plan  will 
report  larger  fuel  consumption  figures  if  this  is  the 
situation . 
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2.  The  model  assumes  constant  wind  factors  over 
each  segment  between  stages,  where  the  computer  flight  plan 
changes  the  wind  factor  as  necessary  at  each  navigation 
aid  or  intersection  between  stages.  This  will  have  a  vary¬ 
ing  effect  on  fuel  consumption. 

In  each  case,  model  fuel  consumption  for  the  air- 
lifter  was  less  than  that  for  the  computer  flight  plan. 

This  is  primarily  due  to  the  computer  flight  plan  carrying 
too  much  fuel  over  the  entire  route.  Reference  to  the  data 
in  Appendix  G  shows  this  to  be  the  case  in  each  of  the 
sample  airlift  scenarios. 

The  question  of  isomorphism — how  true  must  the 
model  be  to  the  actual  system — is  now  addressed.  This 
model  is  not  a  mission  fuel  planning  tool  like  the  computer 
flight  plan;  therefore,  it  is  felt  that  fuel  computations 
need  not  be  as  accurate  as  those  from  the  computer  flight 
plan.  They  should,  however,  be  as  responsive  to  gross 
weight  changes  as  in  the  real  aircraft.  They  also  should 
be  accurate  enough  so  they  do  not  effect  the  optimization 
procedure.  These  example  scenarios,  plus  the  previous 
analyses,  show  that  the  model  meets  these  requirements. 

The  model  is  a  reasonable  representation  of  the  real  system 
and  uses  a  proven  optimization  technique;  therefore,  the 
user  can  feel  confident  that  the  model  is  performing  as 
intended. 
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V  Conclusions  and  Recommendations 


The  primary  objective  of  this  research  was  to 
develop  an  operationally  usable  aerial  refueling  optimiza¬ 
tion  model.  The  various  subobjectives  were  organized  in 
an  effort  to  verify  and  validate  the  results  of  the  model. 
The  overall  conclusion  of  the  analysis  presented  in 
Chapter  IV  is  that  the  AR  optimization  model  found  in 
Appendix  C  performs  as  intended.  Using  dynamic  programming 
techniques,  the  model  determines  optimal  aerial  refueling 
points,  fuel  offloads,  and  tanker  departure  bases,  for  any 
airlifter  mission  scenario.  This  conclusion  is  based  on 
the  following  general  observations: 

1.  Dynamic  programming  is  a  proven  solution  tech¬ 
nique  for  a  problem  of  this  type.  The  Box  optimization  pro¬ 
cedure  was  also  verified  using  a  sample  problem,  the 
results  of  which  are  in  Appendix  F. 

2.  The  results  of  the  experimental  design  and 
analysis  are  reasonable  and  easily  verified.  Response 
variables  varied  as  expected,  from  the  system  analysis 
done  in  Chapter  II,  as  the  control  variables  were  changed. 

3.  Model  output  shows  reasonable  agreement  with 
actual  computer  flight  plan  data.  Fuel  consumption  models 
performed  as  intended.  Also,  model-selected  aerial 


refueling  points  were  surprisingly  close  to  those  actually 
used  on  operational  AR  missions. 

4 .  The  model  shows  a  high  degree  of  isomorphism 
with  the  real  system. 

In  addition  to  verifying  and  validating  the  model, 
the  experimental  design  served  to  answer  some  questions 
on  the  aerial  refueling  problem: 

1.  For  the  scenario  used  in  the  study,  cargo  load 
statistically  had  more  effect  on  total  fuel  consumed  than 
takeoff  fuel. 

2.  There  is  no  general,  direct  relationship 
between  takeoff  fuel  and  total  fuel  consumed.  The  total  fuel 
consumed  is  a  complex  function  of  takeoff  fuel,  tanker  base 
location,  air lifter  route  of  flight,  and  diversion  base 
location.  All  factors  interact  to  determine  the  optimal 

AR  points  and  tanker  offloads.  In  most  cases,  however, 
higher  takeoff  fuel  weights  result  in  lower  total  fuel  con¬ 
sumed  figures  (for  this  scenario). 

3.  For  the  scenario  used  in  the  analysis,  the 
KC-135A  is  more  efficient  than  the  KC-10A  in  the  normal 
rendezvous  mode  of  aerial  refueling. 

4.  The  objective  of  minimizing  total  fuel  con¬ 
sumed  is  not  always  compatible  with  that  of  minimizing  the 
use  of  tankers.  Nevertheless,  the  use  of  the  model  can  be 
highly  beneficial  during  peacetime  operational  mission  plan¬ 
ning  and  contingency  tanker  force  determination. 
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5.  The  location  of  tanker  departure  bases  has  a 
much  greater  effect  on  the  optimal  aerial  refueling  point 
than  originally  expected.  This  fact  alone  may  have  some 
impact  on  future  aerial  refueling  mission  planning. 

Recommendations 

The  following  recommendations  are  submitted  as  a 
result  of  this  research: 

1.  The  planning  of  strategic  airlift  aerial  refuel¬ 
ing  missions  should  include  the  use  of  this  model  to  deter¬ 
mine  optimal  rendezvous  points  and  fuel  offloads.  This 
will  result  in  significant  fuel  savings  to  the  Air  Force, 
and  will  maximize  training  per  flight  hour  for  both  MAC 

and  SAC. 

2.  Model  output  should  be  used  during  the  opera¬ 
tional  planning  phase  of  contingency  planning.  It  can  be 
used  to  determine  tanker  force  requirements  for  a  contin¬ 
gency,  allowing  for  more  efficient  tanker  coordination 
between  MAC  and  SAC. 

Suggestions  for  further  research  in  this  area 
include  the  following: 

1.  To  be  more  beneficial,  the  model  should  some¬ 
how  handle  the  incompatibility  of  minimizing  total  fuel 
consumed  and  the  use  of  tankers.  Reformulation  of  the 
problem  to  consider  both  objectives  would  allow  the  model 
to  be  useful  in  studying  overall  tanker  force  requirements. 


2.  Decreasing  the  fuel  state  step  size  (that  is, 
increasing  NSTEP)  would  result  in  a  slightly  more  accurate 
model.  Also,  increasing  the  number  of  stages,  or  points 
defining  the  airlifter's  route  of  flight,  would  result  in 
output  more  sensitive  to  the  airlifter's  actual  flight 
path.  Making  the  model  computationally  more  efficient  may 
allow  the  user  to  do  both  of  these  without  increasing  com¬ 
puter  execution  time  beyond  a  practical  limit. 

3.  The  "buddy"  rendezvous  method  of  aerial  refuel¬ 
ing  should  be  investigated  and  compared  with  the  results 

of  this  study.  The  results  of  this  study  and  one  on  the 
"buddy"  method  could  be  used  to  come  up  with  an  optimal 
combination  of  airlifters  and  tankers. 

4 .  Wind  factor  data  should  be  used  when  calculating 
tanker  enroute  times.  Furthermore,  a  historical  wind  fac¬ 
tor  data  base  could  be  incorporated  into  the  model,  based 

on  season  and  geographic  location.  This  data  could  be  used 
anywhere  in  the  model  whenever  a  wind  factor  is  required 
and  would  avoid  the  user  having  to  input  this  data. 

5.  The  model  could  be  modified  to  handle  airlifter 
enroute  cells.  This  method  of  airlifter  employment  may 
prove  to  be  the  most  efficient  in  a  large  scale  operation. 

6.  The  model  could  be  modified  to  allow  tankers  to 
depart  with  only  enough  fuel  for  the  mission,  rather  than 
maximum  fuel  load.  This  study  would  result  in  further  fuel 
savings.  Also,  making  the  tanker  fuel  flows  dependent  on 


actual  tanker  gross  weight  rather  than  an  average  gross 
weight  would  further  refine  the  model. 

7.  Finally,  using  actual  computer  flight  plan 
distances  for  both  airlifter  and  tanker  routes  of  flight 
would  increase  model  accuracy. 
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Models 


The  general  model  used  in  deriving  an  aircraft 
fuel  flow  equation  is  the  following: 

a.  b2 

FF  =  bQ(TAS)  (GW)  +  (14) 

where  b  ,  b, ,  and  b„  are  constants,  TAS  is  the  aircraft 
true  airspeed  in  knots,  GW  is  the  gross  weight  in  pounds, 
and  FF  is  the  estimated  fuel  flow  in  gallons  per  hour. 

This  equation  can  be  reduced  to  the  general  log-linear  form: 

ln(FF)  =  ln(bQ)  +b1ln  (TAS)  +  b2  In  (GW)  +  e±  (15) 

from  which  a  multiple  linear  regression  can  be  performed. 
Reference  12  uses  this  procedure  for  determining  a  general 
fuel  consumption  model  for  all  Air  Force  aircraft  collec¬ 
tively;  however,  this  appendix  will  apply  this  procedure 
to  refine  the  equation  for  use  in  the  AR  optimization  model. 

Table  A-l  summarizes  the  data  used  in  the  regres¬ 
sion.  One  model  was  derived  for  both  airlifters  (C-5A  and 
C-141B) ,  and  one  for  each  of  the  tankers  (KC-10A,  KC-135A) . 
Aircraft  gross  weights  were  selected  throughout  the 
operating  range  of  each  aircraft,  and  true  airspeeds  were 
calculated  based  on  the  optimum  altitude  capability  of  each 
aircraft  as  a  function  of  gross  weight.  Therefore,  the 
data  used  assumes  that  each  aircraft  flies  as  close  as 
possible  to  its  optimum  altitude  based  on  its  present  gross 


TABLE  A-l 


DATA  USED  FOR  FUEL  FLOW  REGRESSION 


KC-135A 

GW  (lb) 

Opt  Alt 

TAS  (kt) 

FF (lb/hr) 

FF (gal/hr) 

120,000 

45,000 

470 

7121 

1096 

140,000 

42,000 

470 

8246 

1268 

160,000 

39,000 

4  70 

9216 

1418 

180,000 

36,500 

470 

10444 

1607 

200,000 

34,000 

473 

11537 

1775 

220,000 

32,000 

477 

12892 

1983 

240,000 

29,000 

485 

14265 

2195 

260,000 

26,000 

491 

15839 

2437 

270,000 

24,000 

495 

17067 

2626 

C-141B/C- 

5A 

GW  (lb) 

Opt  Alt 

TAS  (kt) 

FF (lb/hr) 

FF  (gal/hr) 

190,000 

44,800 

433 

7600 

1169 

210,000 

42,700 

433 

9000 

1385 

230,000 

40,800 

433 

10200 

1569 

250,000 

39,200 

433 

11000 

1692 

270,000 

37,500 

433 

11800 

1815 

290,000 

35,900 

433 

12800 

1969 

310,000 

34,300 

436 

13800 

2123 

330,000 

32,800 

440 

14000 

2154 

400,000 

41,500 

443 

14800 

2277 

450,000 

39,500 

443 

16720 

2572 

500,000 

37,000 

443 

18320 

2818 

550,000 

35,000 

445 

20160 

3101 

600,000 

33,000 

448 

22400 

3446 

650,000 

29,000 

4  56 

24800 

3815 

700,000 

27,500 

468 

27120 

4172 

weight.  This  is  what  is  done  in  practice  whenever  pos¬ 
sible.  Also,  the  temperature  deviation  from  standard  was 
assumed  to  be  +10  degrees  Centigrade,  which  again  is  a 
typical  situation.  The  fuel  flow  for  each  data  point  was 
obtained  from  the  applicable  aircraft  performance  manual  in 
pounds  per  hour,  and  converted  to  gallons  per  hour  using  the 
conversion  factor  of  6.5  pounds  per  gallon.  C-5A  fuel  flow 
data  alone  were  used  in  deriving  the  KC-1CA  fuel  flow 
model;  Reference  5  presents  an  analysis  justifying  the  use 
of  C-5A  data  to  estimate  KC-10A  fuel  consumption. 

Regression  Results 

The  Statistical  Package  for  the  Social  Sciences 
(SPSS)  multiple  linear  regression  routine  was  used  in  all 
three  cases,  resulting  in  the  following  fuel  flow  models: 

Airlifters:  ln(FF)  = -6 . 3126+. 83441n (GW) + . 55031n (TAS) , 

or  FF  =  1.81x10_3(GW) ,8345(TAS) *5503  (16) 

KC-10A:  In (FF)  = -11 . 7366+ . 9 79 51n (GW) +1 . 12141n (TAS) , 

or  FF  =  7.995xl0"6 (GW) ' 979 5 (TAS) 1 " 1214  (17) 

KC-135A:  In (FF)  -  -17 .4898+ . 91611n (GW) +2 . 2391n (TAS) , 

or  FF  =  2.54xl0_8 (GW) *9161 (TAS) 2,239  (18) 

The  SPSS  Summary  Table  for  the  airlifter  fuel  flow 
model  is  shown  in  Table  A-2,  and  the  predicted  and  actual 
fuel  flows  are  plotted  in  Figure  A-l.  Statistically,  the 
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TABLE  A- 2 


SPSS  SUMMARY  TABLE— AIRLIFTER 


TAS 

GW 

F  to  Enter  or  Remove 

.13615 

112.97 

Significance 

.719 

.000 

Multiple  R 

.889 

.990 

R  Square 

.791 

.979 

R  Square  Change 

.791 

.189 

Simple  R 

.89 

.99 

Overall  F 

292.5 

Significance 

.000 

model  is  highly  significant,  with  a  coefficient  of  deter- 
2 

mination  (R  )  value  of  .979.  This  indicates  that  98  per¬ 
cent  of  the  variability  in  fuel  flow  is  explained  by  the 
regression  model.  The  high  overall  F  value  of  292.50  would 
cause  a  rejection  of  the  hypothesis  Hq,  where 


31  ~  B2  *  0 

at  least  one  inequality 


confirming  that  true  airspeed  and  gross  weight  together 
contribute  significant  explanatory  power  to  the  regression 
model.  In  a  fuel  consumption  estimating  sense,  Figure  A-l 
shows  reasonable  agreement  between  actual  and  predicted 
fuel  flows . 

Table  A-3  contains  the  Summary  Table  for  the  KC-135A 

model,  and  Figure  A-2  shows  a  comparison  between  actual  and 
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TABLE  A -3 


SPSS  SUMMARY  TABLE— KC-135A 


TAS 

GW 

F  to  Enter  or  Remove 

46.25 

1446.73 

Significance 

.000 

.000 

Multiple  R 

.897 

.999 

R  Square 

.805 

.999 

R  Square  Change 

.805 

.194 

Simple  R 

.897 

.996 

Overall  F 

3712.58 

Significance 

.000 

2 

predicted  fuel  flows.  Here,  R  is  .99,  indicating  that  the 
model  is  able  to  explain  99  percent  of  the  variability  in 
fuel  flow.  The  extremely  high  F  value  of  3712  confirms  the 
result  as  a  highly  significant  model  in  a  statistical 
sense.  Figure  A-2  shows  extremely  good  agreement  between 
actual  and  predicted  fuel  flow  values. 

The  Summary  Table  for  the  KC-10A' model  is  given 
in  Table  A-4 ,  also  indicating  a  highly  significant  model. 
Again,  99  percent  of  the  variability  in  fuel  flow  is 
explained  by  the  model.  Figure  A-3  shows  good  agreement 
once  again  between  predicted  and  actual  values. 


i 
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TABLE  A- 4 


SPSS  SUMMARY  TABLE  — KC-10A 


TAS 

GW 

F  to  Enter  or  Remove 

16.85 

1207.16 

Significance 

.015 

.000 

Multiple  R 

.861 

.999 

R  Square 

.741 

.999 

R  Square  Change 

.741 

.258 

Simple  R 

.861 

.998 

Overall  F 

2339.43 

Significance 

.000 

Validity  of  the  Model 

Yi  =  30  +  +  02X2  +  e±  (19) 

Given  the  above  general  model,  the  linear  regres¬ 
sion  procedure  assumes  the  following  to  be  true: 

1.  It  assumes  that  the  mean  of  Y,  given  X^,  X^, 
is  predicted  by  the  above  multiple  linear  regression  model. 

2.  The  error  term,  e^,  is  a  random  variable, 
normally  distributed,  with  its  mean  equal  to  zero. 

3.  also  has  a  constant  variance. 

All  residuals  were  tested  by  the  SPSS  routine  and  found  to 
be  within  two  standard  deviations  of  the  mean.  Although 
this  test  does  not  confirm  normality  of  the  residuals,  it 


of  Fuel  Consumption  Model 
Data — KC-10A 
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does  indicate  a  lack  of  evidence  against  the  normality 
assumption;  therefore,  it  is  assumed  that  the  second  assump¬ 
tion  is  a  valid  one.  The  constant  variance  assumption  of 
the  error  term,  was  checked  by  plotting  the  residuals 
against  FF,  the  predicted  fuel  flow  values.  These  plots 
are  shown  in  Figures  A-4,  A-5,  and  A-6.  In  each  case,  the 
residual  appears  to  be  independent  of  the  change  in  the 
predicted  fuel  flow;  therefore,  the  constant  variance 
assumption  is  a  valid  one. 

The  variance/covariance  matrix  (as  determined  by 
SPSS)  for  each  of  the  three  fuel  flow  models  indicates 
that  the  coefficients  as  determined  by  the  regression  are 
essentially  independent.  The  small  covariance  between 
variables  further  substantiates  the  validity  of  the  models. 
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This  section  explains  in  detail  the  computer  code 
for  the  AR  optimization  model  contained  in  Appendix  C. 
Figure  B-l  contains  a  logic  diagram  showing  how  the  main 
program  DYNAM  interacts  with  the  various  subroutines.  A 
separate  appendix  (Appendix  D)  is  devoted  to  explaining 
the  Box  "Complex"  algorithm. 

The  program  consists  of  a  main  program  (DYNAM) , 
three  general  subroutines  (COMPLX,  CHECK,  and  CENTR) ,  and 
seven  other  subroutines  designed  specifically  for  the  AR 
optimization  model  (FUNC,  CONST,  SCONST,  RETURS,  TRANS, 
DISXRT,  and  DISBAS) .  The  main  program  coordinates  the  sub¬ 
routines  to  provide  the  stage  optimums  and  performs  the 
recursion  to  recover  the  optimum  path.  Final  returns, 
decision  values,  and  state  values  are  printed  from  DYNAM; 
intermediate  results  can  also  be  printed  from  subroutine 
COMPLX  if  the  user  desires.  The  main  program  and  three 
general  subroutines  were  obtained  from  Reference  16,  with 
some -modification .  All  subroutines  were  specifically 
written  for  the  AR  optimization  problem. 

DYNAM  has  several  user-selectable  options  for 
printing  results  and  formulating  the  problem.  The  param¬ 
eter  IPRINT  selects  the  printing  option  to  be  used  accord¬ 
ing  to  the  following: 


Input : 


Fig.  B-l .  DYNAM  Flow  Diagram 


STOP 
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IPRINT  =  0:  Only  print  final  results. 

IPRINT  =  1:  All  search  values  and  tables  are  printed 

IPRINT  =  2:  State,  return,  and  decision  values  are 
printed  after  each  stage  is  optimized. 

The  parameter  IPROB  selects  the  type  of  problem  to  be 
solved: 

IPROB  =  1:  Initial  value  problem;  given  a  takeoff 
fuel  weight,  the  program  will  optimize  stage  decisions 
without  regard  to  the  overhead  destination  fuel. 

IPROB  =  2:  Final  value  problem;  given  a  required 
overhead  destination  fuel,  the  program  will  optimize  stage 
decisions  and  will  determine  takeoff  fuel  weight. 

IPROB  =  3:  Initial  value/final  value  problem; 
given  both  takeoff  fuel  and  required  overhead  destination 
fuel,  the  program  determines  the  optimal  stage  decisions. 
This  option  is  used  in  this  analysis. 

DYNAM  also  has  the  capability  to  optimize  an  objec¬ 
tive  function  over  more  than  one  decision  at  each  stage; 
for  example,  the  model  could  be  modified  to  determine  the 
optimal  type  of  tanker  and  number  of  tankers  at  each  stage'. 
The  parameter  NDECIS  defines  the  number  of  types  of  deci¬ 
sions  at  each  stage.  To  handle  this  modification,  however, 
dimension  statements  and  common  blocks  would  have  to  be 
changed  (see  Reference  16)  . 

The  following  discussion  describes  each  of  the 
unique  subroutines  in  detail. 


Subroutine  SCON ST 


Subroutine  SCONST  is  called  by  the  main  program 
before  the  recursive  analysis  of  each  stage  to  set  the 
minimum  and  maximum  allowable  fuel  states  allowed  at  each 
stage*  In  each  case,  these  constraints  apply  to  the  air- 
lifter's  fuel  state  prior  to  any  refueling  done  at  that 
stage . 


Minimum  State  Bound.  For  stage  1  (overhead  des¬ 
tination)  ,  the  lower  bound  on  the  state  value  is  determined 
by  the  required  overhead  destination  fuel  figure.  For  all 
other  stages,  the  subroutine  calculates  a  fuel  requirement 
for  the  airlifter  to  divert  to  the  nearest  diversion  base 
from  that  stage  (SCONl) .  It  also  calculates  the  fuel 
required  for  the  airlifter  to  fly  the  next  portion  of  its 
route  of  flight  (SC0N2) .  SCONST  then  selects  the  lowest 
of  these  two  fuel  figures  to  be  the  state  lower  bound  at 
that  stage.  This  insures  that  the  airlifter  will  have 
enough  fuel  to  either  divert  to  a  suitable  alternate  or 
continue  to  the  next  stage  if  refueling  is  either  unsuc¬ 
cessful  or  not  accomplished. 

Maximum  State  Bound .  For  stage  1,  the  upper  bound 
on  the  state  value  is  again  the  required  overhead  destina¬ 
tion  fuel.  This  forces  the  main  program  (DYNAM)  to  come 
as  close  as  possible  to  this  fuel  state  during  the  recur¬ 
sive  calls  to  recover  the  optimal  system  solutions,  thereby 
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discouraging  the  offload  of  excess  fuel.  The  lower  bound 
acts  like  a  hard  constraint,  whereas  the  upper  bound  does 
not.  For  the  remaining  stages,  the  subroutine  selects  the 
state  value  as  constrained  by  the  airlifter's  maximum  fuel 
capacity  or  maximum  allowable  gross  weight,  whichever  is 
more  restrictive. 

The  main  program  then,  using  the  state  value  con¬ 
straints,  divides  the  allowable  state  space  into  100  incre¬ 
ments  (as  defined  by  the  parameter  NSTEP) .  Starting  with 
stage  1  and  the  lowest  possible  state  value,  DYNAM 
coordinates  the  rest  of  the  subroutines  to  determine  an 
optimal  decision  value,  or  number  of  tankers  to  be  off¬ 
loaded,  for  each  of  the  100  possible  state  values.  This  is 
done  to  minimize  the  return  function  for  each  stage,  and 
the  optimal  decision  values  and  returns  for  each  possible 
state  are  stored  for  recall  when  the  system  optimal  solu¬ 
tion  is  recovered  by  DYNAM. 

Subroutine  CONST 

Again,  the  general  problem  is  to 

Minimize  F (r. ,r_, . . . ,r  ) 
l  z  n 

Subject  To  Gi  <  di(Si)  <  Hi  (20) 

The  subroutine  CONST  determines  the  decision  value  con¬ 
straints  as  a  function  of  the  fuel  offload  requirements. 

The  decision  value  for  a  stage  represents  the  number  of 
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tankers  offloaded  at  that  stage.  Therefore,  the  con¬ 
straints  are  defined  in  CONST  as  fractional  parts  of  the 
tanker’s  available  fuel  for  offload,  or 

=  (Fuel  Requirement) /FTNKi  (21) 

where  is  the  lower  constraint,  and  FTNK^  is  the  maximum 
amount  of  fuel  one  tanker  can  offload  at  stage  i. 

Lower  Constraint.  CONST  first  calculates  the  fuel 
required  by  the  airlifter  to  fly  to  the  next  stage  (FI) . 

It  then  calculates  the  fuel  required  to  fly  from  the  next 
stage  to  the  closest  diversion  base  (F2) .  If  the  sum  of 
FI  and  F2  is  greater  than  the  state  value  presently  being 
considered,  the  lower  bound  becomes 

[F1+F2- (STATE  VALUE) ]/FTNK  (22) 

If  not,  the  lower  bound  is  zero.  This  insures  that  if  the 
state  value  presently  under  consideration  will  not  allow 
flight  to  the  next  stage  and  diversion  at  that  stage  (if 
AR  is  unsuccessful) ,  the  airlifter  will  be  forced  to  refuel 
to  at  least  meet  this  requirement.  For  stage  1,  the  lower 
constraint  is  set  to  zero,  since  it  would  be  meaningless 
to  refuel  overhead  destination. 

Upper  Constraint .  The  upper  bound  on  the  decision 
variable  is  also  set  by  either  the  maximum  fuel  capacity, 
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or  the  maximum  gross  weight  of  the  airlifter,  whichever  is 
more  restrictive: 


Hi  =  (XMAX-SLOW)/FTNKi  (23) 

where  XMAX  is  either  the  maximum  fuel  capacity  or  the  maxi¬ 
mum  fuel  weight  as  limited  by  maximum  gross  weight,  and 
SLOW  is  the  lowest  possible  state  value.  For  stage  1,  the 
upper  bound  is  also  set  to  zero. 

If,  for  a  particular  stage,  the  tanker  bases  are 
so  far  away  that  FTNK  is  less  than  zero,  CONST  returns  a 
value  for  and  of  zero.  This  eliminates  any  possi¬ 
bility  of  refueling  at  a  point  which  is  infeasible  and 
forces  the  recursive  analysis  to  look  for  fuel  elsewhere 
in  the  system. 

These  bounds  on  the  decision  values  are  used  by 
COMP LX  (see  Appendix  D)  to  select  an  initial  set  of  feas¬ 
ible  decision  values  for  the  recursive  optimization  pro¬ 
cedure.  Subroutine  CHECK  is  used  by  COMPLX  to  keep  all 
possible  decision  values  during  each  iteration  within  these 
.limits. 

Subroutine  RETURS 

Subroutine  RETURS  calculates  the  total  return  for 
each  stage  as  a  function  of  the  state  value  and  decision 
value  under  consideration.  The  stage  return  is  defined  as 
the  fuel  required  by  the  tanker (s)  to  fly  from  departure 
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base  to  the  stage  and  return,  plus  the  fuel  required  by 
the  airlifter  to  fly  from  that  stage  to  the  next  one  along 
its  route  of  flight.  RETURS  is  called  by  subroutine  FUNC 
each  time  either  the  stage,  state,  or  decision  value 
changes  during  the  iterative  process.  The  return  for  each 
combination  at  each  stage  is  used  by  COMPLX  during  the 
optimization  procedure,  since  the  objective  function  is 
the  sum  of  the  individual  stage  returns: 

n 

Minimize  F (r, ,r, , . . . ,r  )  =  E  r.  (24) 

1  z  n  i=1  i 


Subroutine  TRANS 

Subroutine  TRANS  is  called  by  DYNAM,  COMPLX,  and 
FUNC  any  time  a  state  transformation  is  required  between 
stages.  TRANS  contains  the  transition  function  for  the 
dynamic  programming  procedure.  The  subroutine  begins  with 
the  fuel  state  at  stage  i+1,  subtracts  the  fuel  required 
by  the  airlifter  to  fly  to  stage  i,  and  adds  the  fuel 
onload  at  stage  i+1  to  determine  the  fuel  state  at  stage  i 
The  transition  is  a  function  of  the  state  values  and  deci¬ 
sion  values  under  consideration. 

Subroutine  FUNC 

This  subroutine  formulates  the  objective  function 


for  each  state  value  before  the  optimization  over  all 
decision  values  process  begins.  COMPLX  is  the  primary 


user  of  subroutine  FUNC,  and  FUNC  itself  calls  RETURS  and 
TRANS  while  determining  the  objective  function.  For  each 
stage,  COMPLX  calls  FUNC  with  the  state  value  under  con¬ 
sideration  and  its  associated  decision  values  as  inputs, 
from  which  a  return  is  calculated  for  each  decision  using 
RETURS.  FUNC  then  uses  subroutine  TRANS  and  the  state 
value  to  calculate  the  corresponding  state  of  stage  n-1. 
Then,  FUNC  searches  through  the  table  of  state  values  for 
stage  n-1  until  it  matches  the  calculated  state  value  from 
TRANS  with  one  of  these  stored  state  values.  Since  the 
stored  state  value  has  an  optimal  decision  and  function 
value  associated  with  it,  based  on  the  recursive  analysis 
up  to  that  point,  the  stored  function  value  for  stage  n-1 
is  added  to  the  return  calculated  from  RETURS  to  form  the 
new  function  value.  The  sum  represents  the  objective  func¬ 
tion  for  stage  n  given  the  state  value  under  consideration 
and  is  the  function  to  be  minimized.  This  is  done  each 
time  the  decision  value  is  changed  during  the  optimization 
procedure  in  COMPLX. 

Subroutine  DISBAS 

Subroutine  DISBAS  calculates  the  great  circle 
distance  in  nautical  miles  between  two  points  using  geo¬ 
graphic  coordinates.  The  following  expression  is  used  to 
calculate  the  angular  difference  in  radians: 


9' 


D  =  cos”'L  [sin  (MACLA)  *sin  (XP1) +cos  (MACLA)  *cos  (XP1) 

*cos (abs (XP2-MACL0) ) ]  (Ref  5:168)  (25) 

where 

MACLA  =  Latitude  of  the  stage 

MACLO  =  Longitude  of  the  stage 

XP1  =  Latitude  of  the  tanker  departure  base 
XP2  =  Longitude  of  the  tanker  departure  base 

The  distance  in  nautical  miles  is  then  computed  by  using 
the  following  relationships: 

1.  One  degree  of  latitude  =  60  nautical  miles. 

2.  One  radian  =  57.2958  degrees. 

D1SBAS  computes  the  distances  between  the  stage 
and  all  tanker  departure  bases  included  as  input.  It  then 
selects  the  base  with  the  smallest  distance,  stores  the 
matrix  index  associated  with  that  base,  and  returns  the 
smallest  distance  as  the  variable  XTNK.  XTNK  is  then  used 
by  DYNAM  to  compute  the  maximum  permissible  fuel  offload 
per  tanker  at  each  stage,  FTNK^ . 

DISBAS  is  also  used  to  compute  the  distance  between 
each  stage  and  all  diversion  bases.  It  also  selects  the 
base  with  the  smallest  distance,  stores  the  matrix  index 
for  that  base,  and  returns  that  distance  as  the  variable 
XDIV.  XDIV  is  then  used  by  CONST  and  SCONST  to  calculate 
the  diversion  fuel  required  at  each  stage. 
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Subroutine  DISXRT 


Subroutine  DISXRT  uses  the  same  relationships  pre 
sented  above  to  compute  the  distance  in  nautical  miles 
between  each  point  on  the  airlifter's  route  of  flight. 
DISXRT  returns  the  distance  as  the  variable  XRT 


This  appendix  presents  the  computer  code  for  the 
AR  optimization  model  and  describes  the  variables  used  in 
the  program.  The  general  structure  of  DYNAM  and  its  sub¬ 
routines  were  obtained  from  Reference  16.  DYNAM  has  the 
built-in  capability  to  handle  more  than  one  decision  vari¬ 
able  per  stage.  It  also  has  three  options  for  types  of 
problems  to  be  solved  and  program  output,  all  of  which  are 
described  in  Appendix  B  and  are  selectable  by  the  user. 
These  options  were  left  in  the  program  to  allow  the  user 
some  flexibility  in  applying  the  model  operationally. 

DYNAM  is  written  in  ANSI  Fortran  V.  For  the  20 
stage  model,  the  program  requires  65,000  words  of  core 
memory  on  the  Control  Data  Corporation  (CDC)  6600  computer. 
Execution  time  was  between  120  and  140  seconds  for  the  12 
stage  scenario  used  in  the  analysis  of  Chapter  IV. 

Description  of  Variables 

ALPHA  -  Reflection  factor  (see  Appendix  D) . 

BETA  -  Convergence  parameter  (COMPLX) . 

DECIS  -  Decision  variable. 

DECMAX  -  Optimal  value  of  decision  variable  for  each 
stage . 

DELTA  -  Explicit  constraint  violation  correction. 

DIVLA  -  Latitude  of  diversion  base. 

DIVLO  -  Longitude  of  diversion  base. 
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F 

FMAX 

FTNK 

G 

GAMMA 

H 

I 

IB  AC 

IEV1 

IEV2 

IOPT 

IPRINT 

IPROB 

IT 

ITMAX 

K 

Kl 

KODE 


Value  of  objective  function  (FUNC) . 

Optimal  return  for  the  entire  system. 

Maximum  fuel  offload  capability  of  one  tanker 
at  each  stage. 

Lower  decision  value  constraint  (CONST) . 

Convergence  parameter  (COMPLX) . 

Upper  decision  value  constraint  (CONST) . 

Complex  point  index  (COMPLX) . 

Stage  index  during  backward  stage  recursion 
(I-BACkwards) . 

Index  of  point  with  lowest  function  value 
(COMPLX) . 

Index  of  point  with  highest  function  value 
(COMPLX) . 

Code  for  sign  of  function  at  optimum.  IOPT  =  -1 
for  minimum. 

Code  to  select  printing  of  intermediate  results. 
Code  for  defining  type  of  problem. 

Iteration  number  (COMPLX) . 

Maximum  number  of  iterations  for  stage  optimiza¬ 
tion  (COMPLX) . 

Number  of  points  in  the  Complex. 

Do  loop  limit  (COMPLX) . 

Key  used  to: 

1.  Determine  if  implicit  constraints  are  pro¬ 
vided  (COMPLX) . 

2.  Determine  if  high  or  low  constraint  is 
desired  (DYNAM) . 

3.  Determine  if  system  output  is  desired  (DYNAM). 
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t 


LOAD 

M 

MACLA 

MACLO 

NALFTR 

MBAS 

NDECIS 

NDIV 

NI 

NO 

NSTAGE 

NSTEP 

NTNKR 

OPWT 

OPWTT 

R 

RET 

RETMAX 

S 

SCON 

SHIGH 

SLOW 

SMAX 

SN 

SO 


Air lifter  cargo  load. 

Number  of  constraints  on  decision  variables  at 
each  stage. 

Latitude  of  points  along  airlifter  route  of  flight 
Longitude  of  points  along  airlifter  route  of  flight 
Type  of  airlifter. 

Number  of  possible  tanker  departure  bases. 

Number  of  decision  variables  at  each  stage. 

Number  of  possible  diversion  bases. 

Input  unit  number. 

Output  unit  number. 

Number  of  stages,  or  points  along  the  airlifter* s 
route  of  flight. 

Number  of  intervals  that  the  range  of  state  vari¬ 
ables  is  divided. 

Type  of  tanker. 

Operating  weight  of  the  airlifter. 

Operating  weight  of  the  tanker. 

Random  number  between  0  and  1. 

Value  of  return  variable  during  optimization. 
Optimal  value  of  return  function  for  each  stage. 
Value  of  state  variable  during  optimization. 

Value  of  state  variable  limit  (SCONST) . 

Upper  bound  on  state  value. 

Lower  bound  on  state  value. 

Optimum  state  variable  values. 

System  input  state. 

System  output  state. 
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TAS 


Airlifter  true  airspeed. 

TAST  -  Tanker  true  airspeed. 

TNKLA  -  Latitude  of  tanker  departure  base. 

TNKLO  -  Longitude  of  tanker  departure  base. 

WF  -  Average  wind  factor  between  stages. 

X  -  Value  of  decision  variable  during  optimization 

XC  -  Centroid  (COMPLX--see  Appendix  D) . 

XDIV  -  Great  circle  distance  from  stage  to  closest 

diversion  base  (nautical  miles) . 

XTNK  -  Great  circle  distance  from  stage  to  closest 
tanker  base  (nautical  miles) . 

XTNKM  -  Index  of  closest  tanker  base  from  each  stage. 


PROGRAM  DYNAM  (INPUT  ,  OUT  F  UT ,  T  APES  t  = INPUT ,  T  APECU  =  CUT  PUT) 
HA  IN  LINE  PFOGPAH  FOR  DYNAMIC  PROGRAMMING  ALGORITHM 


THIS  F?OGkAM  DETERMINES  T  HE  OPTIMA;.  AERIAL  REFUELING 
POINTS  AND  FUEL  OFFLOAD  AT  EACH  POINT  TO  MINIMIZE  THE 
FUEL  CONSUME 0  PY  BOTH  AIRLIFT  Er  AND  TANKEF,  IT  ALSO 
DETERMINES  T HF  OPTIMAL  T  A  NKF  *  PASE  FOP  EAC^  PE  FUELING. 
”N CUT  POpMfcT  IS  FRuEFIELD.  THE  PROGRAM  RtDUIRES  THE 
FOLLCWING  DATA  ARRANGED  BY  CAF.U  AS  FOLLOWS 


CARD  #1  -  NUMBER  OF  POINTS  DEFINING  AIRLIFTED  ROUTE 
OF  FLIGHT. 

2  -  CCOR CINATES  ( DEGi  E t S. OEOI MA L)  OF  EACH  POINT, 

BEGINNING  WITH  LEVEL  OFF  POINT  AND  ENDING 
WITH  DESCENT  POINT. 

7  -  NUMBER  OF  TAN KEP  BaSlS  TO  3F  CONSIDERED  (11  MAX). 

«•  -  C  COP  CINATES  OF  TANKER  BASES. 

5  -  NUMBER  OF  AVAILABLE  AI\LI FT  ER  OIVEPSION  BASES  (1.  MAX). 

f>  -  COORDINATES  OF  DIVERSION  BASES.  THIS  SHOULD  INCLUDE 
DFPAFTURE  AND  DESTINATION  IF  THEY  QUALIFY  AS 
DIVERSION  BASES. 

7  -  WIND  FACTOR  FOR  EACH  LES. 

6  -  TAKEOFF  FULL,  RLOU1RLD  OVERHEAD  DESTINATION 

FUEL,  AIRLIFTER  CARGO  LOAD  (ALL  IN  POUNDS). 

9  -  TYPE  OF  AIRLIFTEF 

1  -  C-SA 

2  -  C-l«13 

i'*  -  TYPE  OF  TANKER 

1  -  KC-135A 

2  -  KC-1,;A 

THIS  PROGRAM  WILL  HANDLE  INITIAL  VALUE  (GIVEN  TAKEOFF 
FUEL),  FINAL  VALUE  (GIVEN  REQUIRED  OVERHEAD  DFSTINATION 
CU  E  L)  ,  OR  BOTH  INITIAL  VALUE/FINA.  VALUC  P^OPLtMS. 
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DIHEhSION  X  ( 2'  ,1)  ,  E(2),  CECNAX(2  ,1),  NTCT(2  ),  SMAX  (20  , 

1  0><3),  H(3),  XC  (2),  RETNAX(21>,  DUK9(2L),  DUH91  (20  , 

2TNXLM1),  TNXLO(lw),  OIVLAd.),  DIVLOdO,  XTNK(2!  ) 

00*MCN  N,  S(2f  ,1  1)  ,  IP,  NET  (2(  ,  ivl)  ,  SO,  OECI S  (2C  ,i,  ILi)  , 
lr0P(  2(  »  11)*  NI07N,  NTOr  Ml,  lOPr,  TAS,  OPWT,  LOAD,  XDI  V  ( 2'  )  , 
2FTNK  (20  ,  X  T  T  (  2  :J )  ,  NSTEP,  NSTAGE,  SLOW,  WF(2L),  XMAXFA, 
3XMAX-T,  XMAXGW 

INTFGE9  GAMKA,  X’-JKH(2D,  X01VK2O 
RE  AL  LOAD,  t*  ACL  A  (  ?*' )  ,  MACLO(cl) 

mi  =  r- 

HO  =  fB 

OPTI'iT-'ATILf  P ARA 1ET  EPS  AND  INI T  I  i  tl 7  AT I 0 N.  THESE  PARAMET  EDS 
ARE  THF  ONLY  ONES  CHANGEABLE  BY  THE  USER 

IP  El  NT?  ?:  -  ONLY  PRINTS  FINAL  RESJLTS 

1  -  ALL  SEARCH  VALUES  A  ND“T  A  9L  ES  ARE  PP.lNTfO 

2  -  ONLY  VALUES  AND  TABLFS  ARE  PRINTET  AFTER 

EACH  STAGE  IS  OPTIMIZED 

IPPO’’  J  1  -  INITIAL  VALUE  Pf-OFLEM 

2  -  FINAL  VALUE  FJOBLE* 

3  -  INITIAL  VA.UE/FINA.  VALUE  PROBLEM 

IT  MAX  =  15 
IPdNT  =  !"’ 

IP°Oc  =  3 
TO  F'r  =  -1 
AL  CHA  =  1,3 

beta  =  r . ::  1 
GANNA  =  5 
DELTA  r  i,:  E-t  6 
FMAX  =  '  ,  r 
N  =  i 
*  -  1 

NOFCTS  =  1 
CALL  PANSET  (13) 

DATA  INPUT 

READ  r , NSTA  GE 

DO  l  f  1  =  1, NST  AGE 

RE  *0  *,  KACLA(I),  MACLO(I) 

PR  IN’*-  *,  NACLA(I),  MACLO(I) 

(E  CONTTNUE 

READ  *,  NBASE 

DO  '■'8  1  =  1,  NBASE 

r.EAP  *,  TNKLA(I),  TNKLO(I) 

°r  TNT  TNKLA(I),  TNKLO(I) 

u  »,T  *  fit  J  r 
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RE'D  *,NDIV 

oo  •  o'  :=i,KPTv 

RE  AO  *,  DIVLA(I),  OIVLO(I) 

PRINT  *,  OIVLA(I),  OIVLO(I) 
l  c'ti  CONTINUE 

DO  ft  1=1 ,  NSTAG E-i 

RE  An  *f  MF ( I ) 

PR  ~NT  ' ,WF(I) 
r  •  l  CONTINUE 

REAP  *,SN,  SO,  LOAD 
OR It|T  *  ,  SN,  SO,  LOAD 
PESO  *,KTNKF,  NALFTR 
PR  IN"  n»NTNKR*  NALFTR 

SO  2  =  SO 
IPWGLO  =  1PFINT 
WRIT*.  (NO,  .15) 

'  i  *  FORMAT  (1H1, 1.  X,2  IHOYNAMIC  PRGGRA  W -ILNG  PROCEDURE  > 

Ei  .  WRITE  (NO,-  23)  NSTAGE 

'23  FORMAT  (//,2X, 16HOPTIMIZATION  Or  ,I2,39H  STAGE  INITIAL  AND  FINAL  V 
1ALUE  PROBLEM  ) 

INITIALIZE  AIRCRAFT  PARAMETERS 

IF  (  NALFTR.  EO.l)  THEN 
TAS  s  i'  E  •  ... 

OF  V»T  =  3'Vi  ' 
xmaxfa  =  3ie  I.-- 
X^AX'W  =  7  jl  c  l'. 

EL  SE 

TAS  =  ff32.f. 

OP  WT  =  ll.bgpf; 

XMAXFA  =  1E3  35  2 
XNAXGW  =  3231. 

ENCIF 

IF  (NTNKR.EQ.l)  THEN 
TAST  =  479. <r 
XM'AXCT  =  1  (■  3 1. [  - 
0**WTT  =  i  L  !  '• 

EL  Sc 

TAST  =  <*79.  f 
XM  AXCT  =  3hI 
0°WTT  =  2  A  S  A  7  A 
ENCIF 


USE  SUP  F  CUT  1 NE  S  TO  CALCULATE  DISTANCES  FROM  COORDINATES 
CALL  DTSXRT  (H  ACL  \  ,  MACLO  ,  XFT  ,  NS  T  A  GE) 

OA  LI.  PTSF  AS  (MAfL  AjHACLO,  TIKlA,TN<.0,  NSTAGE,  XT  NX,  N  BASE  ,  XTNKM) 
LL  n!S*V,E  (M  ACs.  > ,  MACL  0 , 01  L*- ,  01  Vl  0,  NOT  M-  1  ,  XD1/,N0IV,  XU  17  '  ) 


THIS  ROUTINE  REVISES  INDEX  ORDER  FDR  RECURSIVE  ANALYSIS 
(STARTING  AT  STAGE  1  -  DESTINATION) 

DO  7'  Isi,NSTAG£-i 
OUMB(I)  =  XFT(I) 

DUMB*  (I)  *  WF(  I ) 
t:  CONTINUE 

DD  ti.  I  =  i,NSTAGE-l 
XRT(I)  =  DUN0(NST AGE-I) 

WF  (I )  =  DUhpi(  NSt  AGE-I) 

M  CONTTNU£ 

no  72  1  =  1  ,NSTA  GE 
DUMB  (I)  =  XT NK  (I) 

OUHPKI)  =  XDIV(I) 

7  2  CONTINUE 

DO  77  I  =  i,NSTAGE 

XTNK(I)  =  DUMB  (NSTAGE-H-I ) 

XOIV(I)  =  OUMc'l  (NST AGE+l-I) 

7  3  CONTINUE 

CALCULATES  AVAILABLE  FUEL  FOR  CrFL  OAD  GIVEN  TALKER  DFFARTUF  E 
BASF.  THIS  IS  DOME  FOP  EACH  FOINT  ON  T  HI  AlRLIFTUv  ROUTE 
OF  FLIGHT, 

IF  (NTNK-.FC.i)  THEN 

HQ  7'*  1*1, NS TAGS 

PRINT  *,'XTNK«,I,»  *  »,XTNKU) 

~TUKtI>  =  XNAXFT  -  ((2,?hE-L  TA5T»  4  2.  2'-  *  C  XUAXFT/2.  ♦  CPWTT)*-*  .Sit  ) 
1 '  (  XT NK(  I)  /T  AST  )  +  C«b  *2)  -  21 L  (  3 
7 A  CONTINUE 
EL  SE 

no  7  7  r  =  l,NSTAGE 

FTNK(I)  =  XNAXFT-(  (/,gg»E-t  ♦  TASr,:^l.i21A»  (XMAXF7/2.+OPWTT) 

1  *  *  •  079  F  )  *  (XTNK(I)  /TAST)  *6.E*2)-U2722 
7  6  CONTt nue 
EN  Cl c 

25 r63  LBS  INCLUDES 

112G3  LBS  -  1  HOUR  ORBIT 

77<T  LBS  -  AOOITICNAL  FUEL  FDR  DUMB,  DESCENT,  APPPOACH 
AND  LANOING, 

6!C'  LSS  -  A 5  HINUTES  HOLDING  if  DESTINATION. 


24  72  2  L  0S  INCLUDES 

1  7722  LES  -  1  HOUR  ORBIT 

1  29" '  L  "S  -  ADDITIONAL  FUEL  FOR  CLIMB,  CESCrNT,  APPROACH 
AN ID  LANDING. 

11E"  LBS  -  4'-;  MINUTES  HOLDING  AT  DESTINATION 
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RECURSIVE  ANALYSIS  OF  STAGES  STARTING  AT  STAGE  1 


1  NS  =  N 

DECIS(N,1,1>  =  1.  1 

<  =  ?*Notc:s 

ilDEC  =  NDEC2S 

IF  (  N.EO.NST  AGE)  THEN 
StCW  =  SN 
GO  TO  7b 

INOir 

KOOE  =  f 

CALI  SCONST  (N ,SLOW, KOOE ) 

7  b  KOCF  =  i 

CALL  FDCNST  (N,SHIGH, KODE) 

PRTN'  r , *  SHI GH  =  * ,  S  HIGH 
PRINT  *,*SLGH  *  SSLOH 

PR  TNT  *  ,  #XF 7  #»  K, *  =  #,XRT(L) 

PRINT  * ,  'XPIV' ,N,  *  =  »,XD1V(M 
PRIN*  « , *FTNK' ,N, •  =  «,rTNKO) 

NS  T  Ea  =  It 

ST  IP  =  (S HIGH- SLOW) /NSTtP 
2  IE  NTOT(N)  =  NSTEP  ♦  1 
NTCTN  =  NTOT(N) 

IF  (N.GE.2)  NTCTM1  *  NTOTU-1) 

OEPFCPM  SEAFCC  AT  STAGE  U  FOE  GIVEN  STATE  VALU^ 

DO  1-C  IP=l,NTOTN 

IF  (TPt.InT  .  1}  >3,  28,  26 
?6  IPRINT  =  i 
28  IH1  =  IP  -  1 

S  (  N,  IP)  =  SLOW  ♦  TMl'STEP 

DO  3’  J= 1 , P0  EC 

3  x(i,  j)  =  orcis  (N,  j,d 

CALL  rOKPLX (NOFC,  N,K , IT M A  X, A l FHA , B ET A , G A HM A, DELI  A , X, F, IT  ,IEV2,NO, 
1  I  PRINT  ,R  ,  G  ,  H,  XC) 


FDP(N,TP>  =  10 FT  ■*  FLIEV2) 
no  •>  J=i,NDEC 

*1  DECISCN, J,1F>  =  X(1EV2,J) 

IF  (TFWOLO  -  1)  11.,,  4 bf  4i> 
k  l  IP  FI  NT  =  IPHOLP 
?•  IF  (IF  -  NTCTN)  l*\ ,b5, ib 
f  V  MR  TT  "  ( r.O,  ’  »  I)  N 

.  -  T  {//,:.  y  ,  ?  ,  u,7  APL‘.  OF  V»  LUES  z0<  ST  LG  ,1D 
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i  f  r 


*  {  l 


11' 


12. 


t?P 


1C? 

1  2  M 

13' 
i  U 


1‘ 

ir 

19' 
1  9' 


WRIT  -  (NC,  *  £>) 

FORMAT  (/  >  11HSTATE  VAL  UT  ,  7  X,  l  4  H  FUNCT  2  ON  VALUE, 11X, 
icHDF~ISIGN  VALUES  , /) 

00  f  '  Z 1  =  1,N10T1 

WRIT:  (KOf.f-K)  SCI,  II),  FDPCf  *11)  ,  OCCIS  (N,  J,II)  ,  J  = 
FORMAT  ( 2  X, 1 PF  14.  *„4X,l?EU.fc  ,3  (  4  X  ,  IP  ti4.  6) ) 

CO  NTTNUE 


IF  (N  -  KSTA6F)  11.,  12- ,  12. 
fl  =  *1  ♦  1 
GO  T0  U 


SOLVING  FOR WAR  PS  FOR  OPTIMUM  PATH 

00  2  II=i,NSTAG: 

IB  AC  =  NSTAGE  ♦  1  -  II 
NT CT II  =  NTCT  (I  nAC  ) 

MO  EC  -  NDECIS 

SM  A X  ( TnAC )  *s  S  (ie\C,NTOTN) 

PSTMAXCIBAC)  =  RETUBAC,  NTCTK) 

00  IIP  K= 1 , NDFC 

OE  CM  AX(  Ir  AC  »  K>  =  OECISdPAC,  x,NTOT^) 

IF  (  ]  BA  C  «  Nt  •  NST  AGE )  GO  10  lb 
FM  AX  s  FOF ( I  BA  C , NTOTN ) 

00  l~r  J=1,KT0TN 

GO  TO  C122t3  2i.  ,122)  ,  IPRDP 

IF  (G(IPAC,J)  -  SO)  IF/-,  13c  i  i.  3. 

IF  ( F  r3  ( I  BP  C  » J  )  -  FMAX)  It.  ,  1?.,  13  . 
FM  A  X  =  FDPd  PAC,  >) 

00  V  -  K= 1 , NCf  C 

OF  CM  AX ( I B  AC  i  K)  =  CEC I S  <  I  B  AC  ,  K,  J) 
SMAX(I=>AC)  =  S  (IB  1C,  J) 

RETM  \X( 1° AC )  =  RET (IB AC, J) 

IF  ( •'P=0F.NE.2)  GO  TO  2 

CONTINUE 

GO  TO  ?!  . 

<0  OF  =  ' 

CALL  TP A MS  (InAC,SMAX,OECMAX,SG,XDOE> 
00  nr  J=1 »  NTOTN 

IF  (G(TC£C,J)  -  SMAXdPAC))  19. ,17', 17 
SM  AX  (TnAC)  =  S ( 1 3  A C , J ) 

PCTMAXdPAG)  =  RET  ( I  PAC,  J) 

00  1  K=  1, NOEC 

0ECMAX(1BAC,K)  =  OECISdBAr,  K,J) 

GO  TO  2i 
CO  NTT  HUE 
CO  NT  T  If- 


t NDEC) 


110 


KODE  =  i 

can  TRASS  (IPAC,  RMAX,DtCMAX,SD,<DDE) 

CO  3  'f'  N=  1,  NST  AGE 

NO  EC  =  MOECIS 
00  2  iC  1 1  =  1  f  NO  EC 
21?  X ( N, II)  =  C  E  CM  AX  (N,II) 

Ip  =1 

CALL  P-TURS  (X  )  N) 

3:-;  CONTINUE 

WRITE  (NO,  ii>  FMAX 

I  1 1  FORMAT  (///,2X,2i-HSYSTLM  OFT  I  HUH  RETURN  =  ,lP£lfe.6> 

WRITE  (SO,  12) 

t 12  FOnM AT  (//, 1 tX , 21 N MAXIMUM  STAGE  RETURNS  ) 

00  A  r  I  =  i,f  Sr/t6E 
fc?1  NR  IT"  (NO,  .  13)  I,  kZTMAX  { I) 

M3  FORMAT  ( /  ,  2  X  ,(  HST  AGE  ,12,1  H  RETURN  =  ,i=>MC.8) 

WRITE  (NO,MA) 

MA  FOFMAT  <//,!: X,17H0PTIMUM  DECISIONS  ) 
no  f-'r  J=1 ,  N  ST  AGE 
WO  EC  =  NDEC1S 

WR IT r  (NO,  it)  J,  (J,  I,  X(J,I),  I = 1, NDEC) 

MC  FOTHAT  (/,2X, ONSTAGE  ,  12  ,  OX,  2  O X ,  ?HX  (  ,  II ,  1H,  ,  13 ,  U  H)  =  ,lFElfc.S)) 
3  1  CONTINUE 

WRITE  (  NC  , '  18)  SO 

Tie  FORMAT  (7,2X,2FHTME  SYSTEM  OUTPUT  3TATr  =  ,1PE16.B) 

WSTMi  =  NST  AGS  -  1 

00  r,  r  j j= i ,  nstmi 

WRITE  (NO,  10)  JJ,  3MAX(JJ) 

MO  FOFHAT  (/,2X,2rHTWE  INPUT  STATE  TO  STAGE  ,I2,3H  =  ,lPE16.8) 

CO  NTT  NUE 

SN  =  SMAX (NSTA  GE) 

WRIT-  ( NO , ■ 1?)  SN 

i  17  FORMAT  (  /,2X,23HTHZ  SYSTEM  INPUT  STATE  =  ,1PE1E.6) 

WRIT-  (NO,  I  !  1) 

!  1  FORMAT  (////, i‘X,38H AERIAL  REFUELING  OPTIMIZATION  PROGRAM  ) 

WF  IT  E  (  NO,r  i  2) 

5‘2  format  </,2-  X,  13H3UWMAPY  OF  FESULTS  ) 

NR  IT  r  ( NO  ,  E  ?  Z) 

c  :  3  FORMAT  (/////,  1GX,27HA1RLIFTC:R  ROUTE  OF  FLIGHT?  ) 

RAD  =  07.  2R!  77951 
TO  r  A  1  =  1,  NS1  AGE 

^L  %  <  I )  =  MAC  LA  ( I  )*  RAO 
MACIC(I)  =  N AC  LO ( I ) *RAD 
:  :  <  CONTINUE 

CO  F  r  1  =  1, (.PASS 
TMKLA  CO  =  T  NX  LA  ( I )  *  F  AD 
TN  XL  ((T>  =  T  NKLO  ( I ) *RAO 
;  OO’.T  "‘''IE 
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03  P  ’*  Is  i,  I-  01  V 

II  VIA  (I)  =  CiVLA  ( t )  *  RAD 

niVLC(I)  =  DIV  L0 ( I ) v  RAD 

-  r-  co  nt  :  u*  i  e 

WRIT"  (Nu,£.  /) 

“7  FO r M AT  <///,l  X,IHSTAGE,:,X,iiHC.33R0INATLSfr  X,11HWIND  FACTOR, i-X, 

1  i  .HOJSTANCr  ( f-M)  ,/) 

XR"(NSTAGE)  = 

WF(NETAGE)  =  ' 

30  p  t  1*1,  NS  TAG? 

IB*  C  =  NSTAGf-  +1-1 

WRIT"  (NO,:  9)  I,  MACLA(IBAC)  ,  M.ft  CL  G  ( 1 3  AC)  ,WF(I)  ,XKT(1) 

'  9  F  0  R  M  A  7  (/,12X,I2,  1X,F6.2,  3X,  F  ?  .  2  ,  E  X  ,  F  5 . 0 , 1  X,FV.2) 

8  CONTINUE 

NR  IT  E  (  K  C  ,  E 1  ) 

I  FORMAT  (//// ,1*  X,  2<  HTAK<£F  DEPARTJRE  BASES*  ) 

WRIT"  ( NO ,£11) 

II  FQ-HAT  (/// , 1 7  X, *  H3A  S  E  N  0  .  X  ,  1  2  H  3  3  3  RDIN  A  T  ES  , /) 

03  £12  I  =  1 »  N  BA  S  E 

WRI’"  (NO, 132)  I,  rN<LA(l )  jTNKLcm 

1?  FORMAT  (/,1BX,  12,  5X,F6.2, 3X,  F7.2) 

12  CO  NT  INUE 

WRIT  l  (NO,:  1*  ) 

1 L  FORMAT  (////,  ltX,?tHAIPLIFTEf  riVERSION  P  A  S  £S  ) 

WR IT  E  (NO,'.  IE  ) 

IE  format  (/// ,  16X,*M8A$E  NO  •  ,fc  X  ,  1 2H  3  D  3  RDI N  A  TE  S  , /) 

DO  £  IB  1  =  1, NOI V 

WRITE  (N0,M7)  1 , OIVl  A  (I)  ,CI  VLO  ( I  > 

17  FORMAT  (/,10X, I2,FX,F6.2, 3X, F7.2) 
it  CONTINUE 

IF  (  SAL  FTR.  £0.  1 )  THEN 
WRITE  (NC,:ib> 

tl  FORMAT  (////, 1EX, 16HTYPE  A IF LIFTERS  ,  16X, SHC-pA  > 

ELSF 

WRIT"  (  NG  ,  >  i  9) 

19  FORM.\t  (////,  ltX,  16HTYPE  AIR  LIFTS?!  ,  16X ,  7  HC  -1  ’•  IB  ) 

FNCIF 

WRIT  (NO, £20  LOAD 

l  FORMAT  (/  ,  if  X,  3  2WCARG0  LoACS  ,2’,  X,  '7,  >:  ,  2X,  r  HLPS.  ) 

WRITE  (NO,:  21)  SW 

21  FORMAT  (/ ,lf X, 1AHTAKE0FF  FUELS  ,13<,f7.  ,2X,EHL°S.  ) 

WR  TT  E  (  HC  ,  v  2  2)  SO  2 

22  FORMAT  (/,lf  X,  32HK0D  OVERHEAD  DESTINATION  FUELS  ,  F7  •  I  »  2X,  £HL  BS  • 
IF  (NTNKR.er.il  THEN 

WRITE  (NO, £23) 

27  FOrMAT  (////, IC.X,  13HTYPE  TeNKt-S  ,  t  9X  ,CHKC-i  3f;A  ) 

EL  SF 
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WRIT'  (NC,t2i.) 

ZU  FO  rM  AT  <////, ltX,13WTYPE  TANKER t  ,i9X,7HK01fA  > 


WRIT!  (NO,«  Zb) 

2*  FORMAT  (/////,  1AX,3HN0.  ,5  X,rW0F“L0A0,  rx,i;  HTOTAL  FUF.l,EX, 

1  9 HCT  VERSION  ,3 X,c  WTANKER  »  **X» L  HF  JEL  ) 

WRI’E  (NO ,5  cl  ) 

2 6  PORuAT  (/V;XfEHSTAGE,3X,7HlAhKE*S,  *X,  5H(LPS)  ,AX,14HC0NSUKED  (L9S) 
1  5  X,-  H9ASEt7  X,t  H3ASE  ,DX,  SHS1  ATE  ,//> 

DO  rT7  1=  i»  NST  AGE 
19  AC  =  NS7  A  GF+i-I 
FLOAT  =  X(1,1)'FTS<(I) 

WR  IT  E  <  NO ,  i;  28)  I,  X  ( I » 1)  ,  FLGA  D  ,  P  •  T  “  A  X(  I)  , XOIV K TRAC)  , XTNKM( I BAC)  , 

1  SHAX(I) 

2  6  FO  CM  \T  </,oX,I2,rX,F^,2,J>X,F7.l,?X,F7..  ,9X,I2,GX,I2,5X,F7.(  ) 

27  CONTINUE 

WRIT-  (N0,r?9)  FMAX 

29  FORMAT  (/////,  25X,29HOPriFUF  TCF  A.  FJEL  CONSUMED*  ,3X,F8.v, 
i  2X.9HLBS.  ) 


3U  -POUTIM  COMPLX  IT  MA  X,  AL  ? HA ,  B ET  A ,  GA MM A ,  DcLT A  ,  X ,  F ,  IT , It  V2 , 

1  NC,  TPRINT»P  ,G,H,XG> 

COORDINATES  SPLC14L  PURPOSE  SUBROUTINES 

ARGUMENT  LIST 

IT  =  ITERATION  INDEX 

IE  VI  =  INDE  X  OF  POINT  WITH  MINIMUM  "UNCTION  VALUE 

I£V2  =  INDEX  CF  POINT  WITH  MAXIMUM  FUNCTION  VALUE 

I  =  POINT  INDEX 

KDCC  =  CONTROL  KEY  USED  TO  PL!  ERMINE  IF  IMPLICIT  CONSTRAINTS 
ARE  PROVIDED 
K1  DO  LOOP  LIMIT 

ALL  CTHEPS  PREVIOUSLY  DEFINED  IN  MAIN  LINE 

niMFKSTON  X(K,M),  R(K,N),  F(K),  SH),  H(M>,  XC(N) 

IN  T EG ER  GAhMA 

IT  =  1 
KOPE  «  r 

IF  M-N)  2  ,  2fi  ,  1 ' 

It  KD  DE  =  i 
£  CONTI  HIE 

DO  L  11=2,  K 
DO  3  J  =  J  ,N 
3  x<  II,  J)  =  '  . 

L'  CONTINUE 

DO  LR  11=2, K 
DO  AS  J J=  1 , N 

GENERATES  RANDOM  NUMB  EPS  FCK  CALCULATING  COMPLEX 
P(TI,JJ)  =  FANFO 
it  L  CONTINUE 

IF  (IPPINT)  46  ,  4T,  46 
!.  L  WRIT"  (NO,  :  1) 
til  FORMAT  (//,  2X,  l'JHPARAMET  ERS  ) 

WR IT  E  ( NO , 1  f  2)  N,  M,  K,  IThAX,  IC,  A  _  PH A ,  BETA,  GAMMA,  DELTA 
li2  RORM4.T  (//,2X,/H(  =  ,I2,3X,4HM  =  ,I2,3X,'+H^  =  ,1 2»  2X  ,6  HIT  ^  AX  =  , 
ilL,2Y,“HIC  =.  ,  12, //2X,tHALPH*.  =  ,  F  3  •  2 , 5X ,  7HBETA  V  ,Fi:.5,3X, 

2  BH  GA  MM*  =  , 12, 3X,  3HDELTA  =  ,Ei2.6) 

WRITE  (NO,  '3) 

:,'3  FORMAT  (//,  ?X>  i'+H  RANDOM  NUMBERS) 

00  f  |=2,K 

NpTie  (ko,  r  A)  (J,  I,  R(J,I),  1=1, R) 

•.  PffMr  (/,2X,?(2HR(,I2,lH,,l£,l  H)  =  ,FE.t,2X)> 

CO:  T 
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CM.Cll.ftTE  COMPLEX  POINTS  AND  CH£C<  A  3  A I N  S  T  CONSTRAINTS 


48  DO  EE  11=2, K 
DO  5  J=1,N 

I  =  II 

CALL  CONST  ( N, M,K, X, G,H, I ) 

X<TI,J)  =  G(J)  +  P  <  1 1 ,  J)  *  (H(  J)-5(JM 
?  CONTINUE 
K1  =  IT 

CALL  CMtrK  (N,M,K,X,G,H,I,t<CD£,XC,OELTA,<i) 
r F  (TI-2)  tl,  f  1,  NO 
M  I F  (  1 PE I NT)  02,  63,  0  2 
:  2  WE IT  E  (NO,  16) 

1  £  FO  EM  AT  (//,  2X,  3  JHCOORDINATIS  OF  INITIAL  COMPLEX) 
10  =  i 

WRITE  (NO,  19)  (10,  J,  X(IO,J),  J=l,N) 

1 9  FORMAT  (/,3<2X,2HX(,  12, 1H  ,  ,1  2  ,*H)  =  ,1P£13.6)) 

i .'J  IF  ( "[PRINT)  Ob,  6  5,  56 

•:e  writ-:  (no,  i9)  in,  j,  x(ii,j),  j=i,N) 

E.f  CONTINUE 
XI  =  K 
DO  .7  T  =  1,K 

CALL  FUNC  <N,M,K,X,F,I) 

7  ■,  CONTINUE 
XQUN7  =  1 
IA  =  0 

FIND  POINT  WITH  LOWEST  FUNCTION  VALUE 

IF  ( IF?1NT)  72,  3  i,  72 
7  2  WRIT"  (NC,i  21) 

121  FORMAT  ( / , 2X , 2 2H7ALUES  OF  THE  FUNCTION  ) 

WRIT’'  (NO,'  22)  (J,  F(J),  J-i,X) 

’  2?  FORMAT  (/,2<2X,2MF(,I2,',H)  =  ,l»E13.s)> 

B-  IEV1  =  1 

DO  1  f  ICN=2,K 

IF  (FdEVi) -F(ICW))  lai-9  1  l  i  90 
?  IF  VI  =  ICM 
l:  CO  NT  Tf;')E 

FINn  POINT  WITH  HIGHEST  FUNCTION  VALUE 

IEV2  =  i 
00  12f  ICM=2,K 

IF  (  F  (I  EV2)  -F(  ICM)  )  llr,  11'  t  j.2 
1  1  IE  V2  =  ICM 
12  CONTINUE 


CH  CCK  CONVERGENCE  CRITERIA 


IP  <P(TEV2)  -  F(TEVl)  -  (ArS  (P:Ta*P(I  EVl)  )  )  )  1'.,  1  3C  ,  i 
13  KO*'NT  =  1 

go  n  i* 

It-  KOUNT  =  KOUNT  +  1 

IP  (  l-CM'T-GAMKA)  1E~  >  2*.'  ,  V  i 


REPLACE  F 01 NT  WITH  LOWEST  FUNCTION  VtLU? 

1-  CALL  CENT  K  < N,  ", K,  I c V i, I , XC » X , «1 ) 

DO  l-*-'*  J  J-l  »  N 

If'  X(TEVl,  JJ)  -  (  i.  +  AL  PHA)  *  (XC  ( JJ)  )  -  ALPHA4  (X(I  E VI , JJ)  ) 
I  ~  IF'/i 

CALL  CHECK  (  h,  N,  K  ,  X,  G  ,H,  I  ,KC  Ct ,  X  0 ,  DEL  T  A  ,  <1 ) 

CALL  FUNC  (N,M,K,X,F,I) 


REPLACE  NEW  POINT  IF  IT  REPEATS  65  LOWEST  FUNCTION  VALUE 


17'  I E  V?  =  i 

00  19f  ICK=2,K 

:f  (c(Tev2)“F(JCM))  19-::, is  ,ifc 
15:  I:V?  S  ICH 
13:  CONTINUE 

IF  ( I  EV2-I E  Vi)  22i,2'0,22' 

2 5 1  no  21."  JJ  =  1,N 

X(  IEV1,  JJ)  =  (X(IEV1,JJ)  ♦  XC(JJ))72.:; 

21'  CONTINUE 
I  ’  ICV1 

CALL  CHECK  (  N,  M,  K ,  X,  G  ,H,  I , KO  C  E,  XC  ,  OElTA ,  K 1) 

CALL  FUNC  (N,H,K,X,F,I> 

22'  CONTINUE 

IF  (IPRINT)  23  ,  228,  23. 

2  3  WRIT"  (NO,;  2  3)  IT 

‘£3  FQ ®M AT  (//2X  ,1"HITEP.ATI0N  NUMBER  ,[5) 

WR IT  Z  (NC,  2  A) 

!  2A  c0  THAT  ( / ,  2X  ,  ?  HCQOR  D JN A  T  £S  CF  OORREOTEO  POINT) 
WRITE  (NC,  19)  (I  E  VI ,  JC ,  XCEVi,  JO),  JC=t,N) 
WRITE  (NO,  21) 

WRITE  (NC,  22)  (I,  F(I),  1=1, K) 

WRIT;  (NO,  ?b) 

•  2r  PQ  'hat  ( /  ,  2X  ,  2‘  HC  OOP  DINA  V  ES  OF  THE  CENTROID) 
WRITE  ( N  0  ,  ‘  2 1 )  ( J  T,  XC(JC),  JC=l,X) 

'20  FORMAT  ( /  ,  2  X  ,  3  (2HX  ( ,  12, 5  H  ,C )  *  ,  l»EHi  *6,  *X)  > 

228  IT  *  IT  ♦  1 

IF  (TT-ITMAX)  b 
21  ■  <E  TU  R  N 


sU  IrfvCUTl  NE  C £NTR  ( N»  H,  K,  1  tVJ  ,  1 , XC ,  X  ,  Kl) 

DI!-FFSIOK  X(K,M),  XC  ( N) 

">0  2  ■  J  =  1  ,  N 
XC  (J)  =  .  • 

DO  1  TL  =  1»  Kl 
1  '  XC (J>  =  XC(J)  ♦  X (IL, J) 

OK  =  Kl 

2..  XC<J)  =  (XC(J)  -X(IEVlfJ)  )  /  (R.  K-l.  1) 

iiETIRN 
in  r 


SUORiUTIffE  CHECK  <N, M,K, X ,G, H,l, <DDF, XC, DE L7A, Kl) 
ARGIHFNT  LIST 

ALL  ARGUMENTS  rEFINED  I N  MAIN  LINE  4VD  CONSX 
OIMFfFION  X(K,K),  G(M>,  HIM),  XCM) 

if*  <r  =  c 

CALL  CONST  (  N,  H ,  K ,  X ,  G ,  H,  1 ) 

CH  SC  k  AGAINST  EXPLICIT  OONSTkAHTS 
00  r  s  J  =  1,N 

IF  (X(I,  J)-G(J))  2  C  »  2  u ,  3  C 
2  X<  T,  J)  =  G  (  J )  +  DELTA 

GO  TO  5  ■' 

7  IF  (MJ)-X<I,J>)  4r,**0,SC 

4  X< I,  J)  =  H( J)  -  DELTA 

7t}  COK'TTMJ E 

IF  ( KODE)  11 .» li • 

CHFCk  AGAINST  THE  IMPLICIT  CONSTRAINTS 

i  NN  =  N  ♦  1 
DO  l‘-f  J=NH,M 

CALL  CONST  (N,M,K,X,G,H,I) 

IF  ()I!,J)-G(J))  3v,7  \,7' 

T  IF  (M(J)-X(1,J))  3.,  1.;,'  »  1-., 

5  IF  Vi  =  I 
KT  =  1 

CALL  CENT R  < N, m, K ,1 E VI , I , XC , X ,Kl) 

DO  «T  IJ=1,N 

Y(T,JVJ)  =  (X(I,JJ)  ♦  XC  (  J  J) )  /  2  •  *' 

2'  continue 

.  CO '  IT  TM' E 

:f  t*-)  :  i  ,  i :  *,  i  ' 

i  FT'i)f 
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SUBROUTINE  FUNC  ( N,M, K,X , F,I ) 


dimension  x<k,  n,  F(K) 

DIMENSION  SFUNC(21),  XFUNC(2  ,1) 

COMMON  NUHST,  SI2MU),  IF,  FEr(2,Ul),  SO,  CLCISt  2t;  ,1,11), 
1F0P(  ?P»i-.  1)  *  NTOTN,  NTOTtti,  ICFT,  rAS,  OPWT,  LCAO,  XOIV(2C), 
2FT‘Jt<(2'),  XbU2i),  N  STEP  t  NST  AGE,  SLOW,  WF(2i>,  XHAXFA,  XMaXFT 
RE?  I  LOAD 

CALL  RFTURS  (  X,  1 ) 

IF  (NUMST  -  1)  1  ,  U»  Z<: 

If  P(I)  =  RET  (  MJMST ,  IF) 

GO  TO  <59 


Z'  NUMT  s  NUMST  -  i 

KOPE  =  r 

FS  "M 1  =  FDP(NUMT,NTOTMl> 

3FUNC  (NUMST)  =  3(NUMST,IP) 

00  3  J  F  = 1 ,  N 

3f  XF  UNO (NU^ST , JF )  =  X(I,JF) 

CA  !  L  TPARS  (  NU NT,  SFUNC , X F Uf'C  ,  SO,  KD Ot) 

IF  (MIMT.EO«i*^ND»S(NUMT  ,  l)*GT*SFJslO(NUMT)  )  THEN 
FS7M1  =1.;  L  *f 
SO  TO  9  "• 

ENOTc 

no  0  IS-  i ,  NTO  TMi 

IF  (S  (NUMT  ,  I S)  -  SFtlNC(NUMT)  )  FI,  Vwt 
Vi  FSTM1  =  FOP  (  NUMT,  I S) 

GO  To  e: 
pi  CONTINUE 

5(  F  (  I)  =  PET  (  NUN  ST,  IP)  +  -STUi 

99  F(  I)  r  I  OPT  *  F  ( I ) 

RETURN 


SUBROUTINE  f  CONST  (N  ,  SCON  ,  *00  E> 

COMMON  NUHST,  S(2*,i;l>,  IF,  RET ( ? .  » 1* 1) »  SO,  CECIS(2C,i,l  1), 

ifo e(?f,i  i),  ntotr,  NTorni,  icf r,  fas,  ofwt,  lcad,  xdiv<2d, 

2FTNK  ( 2’ ) ,  XFT(n,  NSTEP,  NS7AGE,  SLOW,  WF(2‘j),  XMAXFA , 

3XH  AX  FT ,  XMAXGW 
REAL  LOAD 

IF  ( KOPc)  ir  , ,  1),  20'1 

1C  IF  ( N  *rO . 1)  THEN 
SCON  =  TO 
GO  TO  99 

>:n  ci  f 

SCONi  =  (i.f  iE-3'-7AS**.55i,,3*  (XNAX'A/2  .♦LOAO+OPWT) ** .83AS* 

1  XCTVCN)/  (TAS*WF(N))  )*f  .  E- 

SCON?  =  (l.£iF-3'-TAS'*.5S  3*  (XKAXFA/2.+LCAD  +  OPWT)**  .63A5* 

1  XRT  (M-1)/(TAS+WF(N)  )  )*•>.£ 

IF  (SCON2.L7 .SCONI)  THEN 
SCON  =  SC  ON  2 
EL  SC 

SCON  =  SCONI 
SN  01 r 
GO  TO  99 

2»o  SCON  =  XHA.XFA 

IF  (  (X^AXFA+OPWT  +  LOAO).GT.XMAXGW)  then 
SCON  =  XMAXGW-CPWT-LOAO 
IN  01  F 

IF  (N.FO.l)  THEN 
SCON  =  SO 

£NOir 


59  RETURN 
ENO 
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SUPP^UTIKE  CONST  (N,M,K,X,G,H,I> 

od-msion  x(k,m),  Gm,  h {»-> 

00**MCN  fUMST,  S(2',i'i>»  IF  »  »;r  T  <  2  ’ ,  lv.1)  ,  SO 
lFDF(?r,  i  i)  ,  N70T  ),  NT0TM1,  luPT,  \US,  0»W”, 
2PTNK(2'),  xncg"),  NSTb»,  NbT  AGE,  SLOW,  WF<2 
3XMAXFT,  XMAXGW 
r.EAL  LOAD 

IF  (FTNKCMff'ST)  .LT»l)  TH:  N 
MM)  =  r 
5( 1)  =  ' 

RETURN 

■:n  oi  f 

IF  (MHST.EO.l)  TH£N 

g<  i)  : 

GO  TO  OS 
SNCIF 

FI  =  (1  •Me-?1  7A3-‘  *  •!/!?■'  3  MC.F  W1 +L.  0  M  ♦  X  M  Ax  r  A /2 
1  *  XPT  (MUMS!  -1)  /(  ra>S*WF<NUKSl  -1))  5  *  5.  v 
r  2  =  (1  .?  U-SMAS  3*  (CP  Ml  *„3S  D*XF.AXrA/2 

i  * XOTV(NUMST-l )/(TAS  +  WF(NUPL i -i)  ) )  S.E 
IF  <  (F1+F2)  ,GT.SMUrlST,IP>)  THr  i 
5(1)  =  (F1  +  F2-S (  JUHST  ,1?)  )/F7 NKlNJISr ) 

■rL  Sr 

g ( i j  =  c.r 

-NCIr 

99  M (  1)  =  (XMkXFA-SwOW) /FTNK(tUHST) 

IF  (  ( yM AXFA  +  OF WT  +  LC'A D) « 3T  »XH;.XGW)  THEN 

H(  1)  =  ((XM{.XGW-OPWT-LOS,0)-SIOW)/FTNK(NU'1ST> 

INFI" 


,  oecis(2f»,i,ui>» 

LCAO,  XDI V  (Zj ) t 
),  XMAXFA, 


)  1.8  3  kl> 
)•  ',8 345 


RETURN 

in  o 


SUBROUTINE-  TRANS  n.BAC»SMAX»D£C'4Atf»S0»K00E) 


DI^ctgn  SRAX  (  2 1 )  *  DECMAX(2.,i) 

CO  WON  N,  S  (  20 , 1  1)  ,  IP,  RET  (  2  ,i:i>,  SOI,  DECISC2- ,  i,  l?  1)  , 
lFDe(V,  V  1)  ,  NTOT1,  NTOTM1,  ICFT,  TAS,  OPWT,  LOAD,  XDIV(2t), 
2FTNK<2  >,  XFTC2T),  NSTEP,  KSTAGE,  SLOW,  M F ( 2 i  ) 

TF AL  LOAD 

IF  (KOPE)  11,  1J,  2f 

1  IB  AC C1  =  IRAC  +  1 

SMAX(r*AC)  =  SNAXdlACPl  )  -( (1  »81E-3*TAS**  •  ?5 3’  (DECHAX  (lEACPlyl  > 

1  “  FTNK(IFaCFI)  +S‘MX(IBA:Px)+UFWrfL3AD)'  “,E2Af)  XRT(IBAC)/ 

2  (  TAS+WFdPAC)  )>*  3  #5  0+DE  C  NAX ( IEA3P1, 1 ) *  FTKK ( TBACPi) 

2  SO  =  SMAX(l) 

ES  RETURN 
END 


SU  Fp  TUT  IKE  fETURS  (X,I) 

OIHEASICN  X(  2f  ,1) 

CO  WHOM  N,  S(2>',ii),  IP,  RET  (  21, 1"  1)  ,  SO,  DECI  S(  2l ,  i,  1  U )  , 

1FO  F( 2 • , 1*  1)  ,  NT  OT  N ,  NTOTMl,  10FT,  TAS,  0=>W7,  LOAD,  XDIV(2D, 
2CTNK(  2)  ,  XF  7  (  2  J)  ,  NSTE3  ,  KST  AGE  ,  SLOW,  WF(2H,  XMAXFA  ,  XMAXFT 
-TEAL  LOAD 

IF  (N.NE.i)  THEN 

RET  (N, TP)  =  <I>i7(<(I,i)  +  .9E3'5>*  (  X'nXFT -FT  NK(  N)  )  )♦  ( ( 1,  8iE-3*T  AS 

1  *  2'  (DFWT  +LO ID  *S  <N,iP)+X(I,l)  *■  FTNK(N))V*»83A5>)*XR7  ( N-i )  / 

2  (  TAS+WF  (N-i ) )  )■*!>•!> 

EL  SF 

RET  (  N,TP)  =  IN7(X  (I,1>+.999S  >*  <X«  AXFT -FT'JK  (N  >  > 

ENPT  c 

RETURN 

INC 


hi  t J 
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SU  nROUTINE  DISPAS  (MACLA,  MACL L ,  x:> i ,  XP2  ,  NST  A GE  ,  X2 , NBAS  , X 3 > 

CO MP'JTES  GPEAT  CIRCLE  DISTANCE  -  11NKEF  AND  DIVERSION  BASES 

DI MF  AST  ON  l^FKi  '),  XP2<1(>,  Xl(2i,lj),  X2(2f) 

INTEGER  X 3 ( £ V ) 

REAL  HA  CL  A  (  2  >  ,  M1CL0(2') 

RAD  =  e7. 29M7P51 
DC  1  1=1,  NBAS 

VP  1  ( I )  =  XP1  (I) /RAD 
XP2(I)  =  XP  2(1) /RIO 
1'  CD  NT  XN’UE 

DO  2  1  =  1, NS TAGE 

DO  7  J=1 ,NBAS 

XI  (T  ,  J)  =  Z‘  37.7v$77*ACDS<SIN(MAC.UI))*SlN(XPi(J>>  + 

1  DOS  (MACLA  (I )  )  '  CDTCXPi(J)  )’CoS(A2S(XP2(J)-MACLO(I>))) 

IP  (  I.EO.i)  THEN 
X2  (1 >  =  X1(I,  J) 

X3  (I )  =  J 
GO  TO  7- 
EN  DI r 

TF  (XiCI, J) ,LT .X2(IJ)  THEN 
V2  (I)  =  X1(I,J) 

XT (I )  =  J 
“N  Cl  c 

'  CONTINUE 
!.  CONTINUE 
RETURN 
END 


SUBROUTINE  DISXRT  (HACLA  ,  HALLO, XU,  NSTAGE) 

CO  f*P  UTrS  GREAT  CIRCLE  DI  STANCE/ MR  LIFTER  ROUTE  (NO.  HEMISPHERE) 

PEAL  MA  C  L  A  ( 2  i  )  ,  HACLO(2’),  XP.T(2:> 

RA  C  =  e7 . 2Gr  77  PT1 
00  1  1=1, NSTAGE 

HA  CL  A (I )  =  F  ACLA ( I ) /RAO 
lAr  LD  <I)  =  M  AC  LO ( I )  /KAD 
1  CO  NT  I  Ml  £ 

L  =  NST  AGE  -  1 
00  2  1  =  1, L 

XRT  CD  =  3M7.74G77-AC0S  (SJN  (NACLMI)  )*SIN(MACI  A(I+1))  + 

1  CCG(M*CLA(I  ))~C03(MACLA  (  2+1 ) )  "  0  DS  ( A9S  ( M6  CL  0  (I  +1) -HA  CL  0  (  J )  )  ) ) 

£  CONTINUE 
RETURN 
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The  constrained  Simplex  (Complex)  method,  as  modi¬ 
fied  for  the  AR  optimization  problem,  searches  for  the  mini 
mum  value  of  the  return  function  r(X^,S^)  subject  to  con¬ 
straints  of  the  form: 

G.  <  X.  <  H.  (26) 

l—i—i 

where  the  subscript  i  refers  to  the  stage  number  presently 
being  minimized.  Only  explicit  constraints  are  applicable 
for  the  AR  optimization  program;  however,  the  algorithm 
also  handles  implicit  constraints  on  the  decision  variables 
if  desired.  The  Complex  program  actually  maximizes  the 
negative  of  the  return  function  in  searching  for  a  minimum. 
It  also  assumes  that  an  initial  point  X^,...,Xn  is  avail¬ 
able  which  satisfies  all  stage  constraints;  for  this  appli¬ 
cation,  a  decision  value  of  zero  is  used  for  the  initial 
point . 

The  algorithm  proceeds  as  follows: 

1.  An  original  "complex"  of  K=3  points  is 
generated,  consisting  of  the  initial  feasible  starting 
point  and  K-l  additional  points.  These  points  are  obtained 
by  using  a  pseudo-random  number  generator  and  the  ranges 
of  the  independent  decision  variable,  such  that 


where  R^  is  a  pseudo-random  deviate  uniformly  distributed 
over  the  interval  (0,1). 

2.  The  selected  points  must  satisfy  both  explicit 
and  implicit  constraints.  If  a  point  does  not  satisfy  the 
explicit  constraints,  it  is  moved  to  a  value  of  6  =  1x10  ® 
inside  the  appropriate  boundary.  If  an  implicit  constraint 
is  violated,  the  point  is  moved  half-way  towards  the 
centroid  of  the  other  two  points,  such  that 


x.  •  =  (X.  .(old)  +  X.  W2  (28) 

-L/J  -L  /  J  l/V' 

where  the  coordinates  of  the  centroid  of  the  remaining 
points,  X.  ,  are  defined  by 


i,c 


J- 


(29) 


This  is  done  until  the  implicit  constraints  are  satisfied. 

3.  The  objective  function  is  evaluated  at  each 
point.  The  point  with  the  lowest  function  value  is  moved 
to  a  distance  a  =  1.3  times  as  far  as  the  centroid  of  the 
remaining  two  points  as  the  distance  of  the  rejected  point. 
The  new  position  is  on  the  line  joining  the  rejected  point 
and  the  centroid,  thus: 


X.  .(new) 
1 1  3 


(X.  -  X.  .  (old) )  +  X. 

ifC  i r  j  j_,c  • 


(30) 


(Ref  16:371) 
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4.  If  a  point  repeats  as  the  lowest  function 
value,  it  is  moved  one-half  the  distance  to  the  centroid  of 
the  remaining  points.  If  at  any  time  the  new  point  vio¬ 
lates  one  of  the  explicit  or  implicit  constraints,  its 
position  is  adjusted  as  discussed  earlier. 

5.  The  process  continues  until  the  objective 
function  values  for  all  three  points  remain  within  3  =  .01 
for  y  =  5  consecutive  iterations.  An  iteration  is  defined 
as  the  calculations  required  to  select  a  new  point  which 
satisfies  the  constraints  and  does  not  repeat  in  yielding 
the  lowest  function  value  (Ref  16:368).  The  maximum  num¬ 
ber  of  iterations  allowed  during  each  optimization  is 
defined  by  the  parameter  ITMAX. 

Figure  D-l  contains  a  logic  flow  diagram  illus¬ 
trating  the  above  procedure,  as  adapted  from  Reference  16. 

Subroutine  COMPLX 

Subroutine  COMPLX  can  be  thought  of  as  the  main 
line  program  for  the  optimization  portion  of  the  AR  optimi¬ 
zation  model.  It  coordinates  the  special  purpose  and 
general  subroutines  plus  prints  intermediate  optimal  solu¬ 
tions  if  the  user  desires.  Subroutines  CHECK  and  CENTR  are 
the  general  subroutines  used  exclusively  with  COMPLX  during 
the  optimization  process;  FUNC  and  CONST  are  special  pur¬ 
pose  subroutines  which  are  also  used  by  DYNAM  and  other 
subroutines  (see  Appendix  B) . 
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Fig.  D-l.  Box  ("COMPLEX"  ALGORITHM)  Logic  Diagram 
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The  model  parameters  used  in  COMPLX  are  as  recom¬ 
mended  by  Box  in  Reference  6.  These  parameters  and  their 


values  are: 

1 .  a (ALPHA)  =1.3 

2.  3 (BETA)  =  .01 

3.  ITMAX  =  150 

4 .  Y  (GAMMA)  =  5 

5.  K  =  3 

6 .  6  (DELTA)  =  1x10 ~6 
Subroutine  CHECK 

CHECK  uses  subroutine  CONST  to  insure  all  points 
in  the  "complex"  meet  all  explicit  and  implicit  constraints. 
If  a  constraint  is  violated,  CHECK  applies  a  correction  to 
bring  the  point  back  within  the  constraint,  as  previously 
discussed. 

Subroutine  CENTR 

Subroutine  CENTR  is  used  exclusively  to  calculate 
the  centroid  of  the  vertices.  It  is  called  by  COMPLX  and 
CHECK  whenever  a  centroid  is  required. 
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SYSTEM  RLTUhK 


2.7e7re  2 e* 


STAGE 
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A 

J. 
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P':  TUPS' 

S  T  AG 

!  . 
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? 

FJ  TUP  N 

- 

?  »  2E  <?  ‘-44?^  +  :  4 

STAGE 

3 

F.nup" 

= 

i,7  3-3E)"2^r+.-f 

STAGE 

4 

f\F  TUP  N 

= 

F  r  7Pi  .  PA  FP-+  4 

STAGE 

r 

Rf TUT  u 

r 
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STAGE 

e 

RfTUP  -i 

s 

?  ,  •?  C-3E+  4 

STAGE 

7 

RE  TUP  ’J 

r 

?  t  ■'  9  7  rr'  -8r  +  h 
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e 
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- 

2. 7 f  ,r+  4 

STAGE 

c 

F  tTUF  N 

- 

?.T9  ;"<2  i, 

Fig.  E-l.  Sample  Raw  Data  Output 
(Dover  AFB--Moskow) 


op'iinu"  decisions 


THE  INPUT  STM  r  TO  f  T  t  GE  1 

THE  INPUT  STM!  TC  El  AGE  2 

THE  INPUT  STMT  TO  ST  fc  GE  ? 

THE  INPUT  S7MC  TO  STAGE  >' 

THE  INPUT  STM.  TC  ST,Gr  r 

THE  INPUT  STM  •*  TC  STAGE  « 

THE  INPUT  STtT-.  TO  STAGE  7 
THE  INPUT  ST  ATT  TO  SlfGC  8 


V.C0  3:  -•  31  *E+14 
1.  16tt  2~*E+f  5 
■  * i»  eiiT- ?E+; h 
E.t  ?53<;8V»tn 
21M- E  2  •  lE+f  O 
1.1J72E472E+1 5 
1.3'  b9  i  1 T-3E  +  !.» 
1. El  7b 27  Jf+fg 


THE  SYST  M  INPUT  ST*  T  ;  =  1,?'  .-'  t  E  +r'v 

Fig.  E-l — Continued 


AEllfcl  4F  FUELING  0  CT  T”!  ?i Ti  ON  PROGRAM 


SU^'fP.Y  Or  RESULTS 


AIRLIFT  L' f '  ..OUT-!  OF  FLIGHT  t 


STAG' 

COO  r  01 

•If  'rF~ 

WIND  FACTOR 

d;s*.ai«;e  < 

1 

5:  ,£? 

-  37  .  '‘2 

.  31. 

Cft.M 

2 

FT. 
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A?. 

0  3,03 

3 

5*  .4* 

-2,3? 
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031.13 

*T 

Si 

1?  •  r " 
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617.24 

l 

c7,:.; 

3fi  4  '..i 

59. 

t33.u:, 

6 

0^.3, 

6  . 

013.02 

7 

*.*  *>  ’ 

-O.r  • 
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8 
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77, 

LC+. 73 

<3 

t  ,0' 

72.  ^ 

.  « 

(  .TO 

Fig.  E-2. 

Sample 

Formated  Data 

Output 

(Dover 

AFB- -Moscow) 

rM.'Kt?  DEPAF»TMF:G  pases: 


»ASF  S  3.  C  00r 01  Nf  T ES 


2 


??.?•*  - ,  7  3 

4  7,'.?  ft.  36 


AT?  LIFT -F  31  VC IRTON  ^AStS 


3ft  SS  NO.  C  OOF.  PIN  A  IF  S 


4  .r  . 
*>',2  8 
FI. 23 

TYP€  AI  -  LIrT  :.P.  i 
CARGO  L  3 ► D  t 
"  A  Kl  OFF  FIJ  .  Lt 


7  2  .  •'  5 


C-i>A 
174  .  ;  . 

IS  4 fcU  .  L3S 


Or  OV£‘H-JAC  OFCTINATION  F'Jr.L:  7:R‘i,  _nS* 


TYFf.  TM:K"‘M 


Fig.  E-2 — Continued 
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UJ 


Problem 


This  appendix  contains  a  sample  problem  used  to 
validate  the  general  DYNAM  program  before  it  was  modified 
for  the  AR  optimization  model.  The  problem  was  taken 
from  Reference  14  and  formulated  in  Reference  16.  It  was 
solved  on  the  CDC  6600  computer,  using  seven  seconds  of 

* 

CPU  time . 

The  problem: 


Maximize  F 


3 

Z  r. 

i 


Return  Function 
Transition  Function 
Constraint 

State  Variable  Constraints 
Starting  Point 
System  Input 
System  Output 
Parameters 


r^  =  5.0d^~-id^ 

si  -  V-4di 

0  <  d.  <  S. 

—  1  —  1 

0  <  S.  <  5 
di  =  1.0 
Sn  =  5.0 
SQ  =  0.0 
N STAGE  =  3 
IP ROB  =  1 
ALPHA  =1.3 
BETA  =  .01 
GAMMA  =  5 

DELTA  =  1x10 ~6 

K  =  3 

IPRINT  =  0 
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ITMAX  =  500 


IOPT  =  +1 
NSTEP  =  100 

Answers:  F  =  11.43 

dx  =  2.447 
d2  =  1.190 
d3  =  .747 


Figure  F-l  contains  the  final  output  for  this  problem 
Algorithm  answers  were  totally  correct. 


DYNAMIC  PROGRAMS  I  MG  PROCLDURE 

OPTIMIZATION  OR  3  STAGE  INITIAL  VALUE  PROBLEM 


SYSTEM  OPTIMUM  RETURN  =  i  •  i 4  ' ». •'  963E*f  1 

^  r- 

MAXIMUM  STAGE  RETURN 


STAGE 

i 

RETURN 

i.  t 

w 

cr 

• 

w 

<1 

STAGE 

2 

RETURN 

=  3.1,12 

STAGE 

3 

RETURN 

=  2.  i  8?  079G2E+1 

OPTIMUM 

OECISION5 

STAGE 

i 

X(i,l)  = 

2. 41 741  3C-7E  +  J  i 

STAGE 

2 

X<2,1)  = 

1.204-  E797E  +  3C* 

STAGE 

3 

* 

X(3,l)  = 

8.074337 L  t E-Il 

THE  SYSTEM  OUTPUT  STATE  =  7, 233m  3?  13t  +  •;  3 

THE  INPUT  STATE  TO  STAGE  1  =  »•  .2f  :  t  (  i6Jt*VG 

THF.  INPUT  STATE  TO  STAGE  2  =  :  » 7  (  .  l  .  J  t>  JE+1,0 

THE  SYSTEM  INPUT  STATE  =  3 , *  l  ■  U  l - E+J C 


Fig.  F-l .  Sample  Problem  Output 
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Experimental  Design  Output  Data 


The  output  in  this  appendix  is  organized  into  the 
same  four  parts  used  in  the  Analysis,  Chapter  IV.  Parts 
One  and  Two  contain  the  data  used  in  the  main  experimental 
design  and  includes  an  analysis  of  variance  (ANOVA) ;  Part 
Three  addresses  the  sensitivity  analysis  of  model  param¬ 
eters,  and  Part  Four  contains  the  data  used  in  verifying 
the  model  against  actual  computer  flight  plans. 

Parts  One  and  Two 

Tables  G-l  through  G-3  contain  the  output  data 
from  the  main  experimental  design  found  in  Table  5. 

Included  are  the  minimum  total  fuel  consumed  figures, 
total  offload,  and  the  number  of  tankers  used  for  each  fac¬ 
tor  level.  The  fuel  figures  are  plotted  against  takeoff 
fuel  and  cargo  load  in  Figures  9-14  in  Chapter  IV. 

Analysis  of  Variance  (ANOVA) 

An  analysis  of  variance  (ANOVA)  was  performed 
using  the  Statistical  Package  for  the  Social  Sciences 
(SPSS)  ANOVA  routine.  Only  the  C-5A  total  fuel  consumed 
data  were  used  in  the  analysis  (to  insure  two  observations 
per  cell)  with  independent  variables  TOFUEL  (airlifter 
takeoff  fuel)  and  LOAD  (airlifter  cargo  load) .  The  objec¬ 
tive  was  to  determine  which  of  these  two  factors  had  the 
most  effect  on  total  fuel  consumed.  Figure  G-l  shows  the 
resulting  SPSS  ANOVA  table. 


TABLE  G-l 


EXPERIMENTAL  DESIGN  OUTPUT  DATA 
C-5A/KC-135A 


Load 

(1000  lb) 

T.O.  Fuel 
(1000  lb) 

Fuel 

Offload 

Fuel 

Consumed 

# 

Tankers 

100.0 

46.0 

310565 

474690 

4 

71.0 

278417 

424520 

3 

96.0 

255223 

426395 

3 

121.0 

219071 

423080 

2 

146.0 

205584 

412969 

2 

171.0 

183468 

404135 

2 

196.0 

157403 

398837 

2 

125.0 

46.0 

309275 

448497 

3 

71.0 

290502 

436669 

3 

96.0 

269400 

439602 

3 

121.0 

243354 

438611 

3 

146.0 

222904 

441634 

3 

171.0 

197904 

441634 

3 

196.0 

172904 

441634 

3 

150.0 

46.0 

331030 

451123 

3 

71.0 

309174 

454324 

3 

96.0 

285666 

454  589 

3 

121.0 

255035 

449919 

3 

146.0 

233174 

454324 

3 

171.0 

206343 

424239 

2 

196.0 

184778 

425724 

2 
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TABLE  G-2 


EXPERIMENTAL  DESIGN  OUTPUT  DATA 
C-5A/KC-10A 


100.0 


125.0 


150.0 


46.0 

71.0 

96.0 

121.0 

146.0 

171.0 

196.0 

46.0 

71.0 

96.0 

121.0 

146.0 

171.0 

196.0 

46.0 

71.0 

96.0 

121.0 

146.0 

171.0 

196.0 


331206 

305407 

280407 

255407 

230407 

205407 

181206 

319385 

295312 

270959 

252532 

238364 

207599 

183189 

361522 

336522 

311522 

286522 

261522 

230956 

202146 


512939 

512797 

512797 

512797 

512797 

512797 

512939 

540821 

514078 

467583 

474171 

478258 

525368 

528185 

626615 

626615 

626615 

626615 

626615 

555899 

539997 


TABLE  G-3 


EXPERIMENTAL  DESIGN  OUTPUT  DATA 
C-141B/KC-135A 


Load 

{1000  lb) 

T.O.  Fuel 
(1000  lb) 

Fuel 

Offload 

Fuel 

Consumed 

# 

Tankers 

20.0 

30.0 

206594 

247578 

2 

45.0 

191594 

24  7578 

2 

60.0 

176594 

247578 

2 

75.0 

109869 

234235 

2 

90.0 

146594 

247578 

2 

105.0 

131594 

247578 

2 

120.0 

116594 

247578 

2 

35.0 

30.0 

170660 

289035 

3 

45.0 

155660 

289035 

3 

60.0 

140660 

289035 

3 

75.0 

125660 

289035 

3 

90.0 

110660 

289035 

3 

105.0 

95660 

289035 

3 

120.0 

80660 

289035 

3 

50.0 

30.0 

165125 

316393 

4 

45.0 

156953 

290356 

3 

60.0 

131908 

286195 

3 

75.0 

116817 

253858 

2 

90.0 

107432 

257652 

2 

105.0 

90744 

286581 

3 

120.0 

82502 

261625 

2 

Sum  of 

Mean 

Signif 

Source  of  Variation 

Squares 

DF 

Square 

of  F 

Main  Effects 

. 32184E+11 

8 

.402E+10 

.618 

.752 

TORJEL 

.47147E+10 

6 

. 785E+09 

.121 

.993 

LOAD 

. 27469E+11 

2 

.137E+11 

2.110 

.146 

2-Way  Interactions 

•84180E+10 

12 

. 701E+09 

.108 

.999 

TOFUEL  LOAD 

•84180E+10 

12 

.701E+09 

.108 

.999 

Explained 

•40602E+11 

20 

.203E+10 

.312 

.994 

Residual 

.13667E+12 

21 

•650E+10 

Total 

.17727E+12 

41 

•432E+10 

Fig.  G-l.  SPSS  ANOVA  Table 


The  analysis  of  variance  assumes  the  following 


model : 


where 


Yii  =  y 


ou  + 


&•  + 
3 


(aB)ij  + 


e .  . 

13 


(31) 


calculated  total  fuel  consumed  from  model 
output 


V  =  overall  mean  of  the  total  fuel  consumption 
figures 

a.  =  the  effect  of  cargo  load  on  total  fuel  con¬ 
sumed 


the  effect  of  takeoff  fuel  on  total  fuel 
con sumed 


(cxf5) 


ij 


the  effect  of  interaction  between  the  two  fac¬ 
tors  on  total  fuel  consumed 


e 


ij 


error  term 
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Using  a  significance  level  of  .95,  the  SPSS  output  indi¬ 
cates  that  none  of  the  factors,  including  interaction,  are 
statistically  significant  in  explaining  the  model.  However, 
cargo  load  does  have  a  much  greater  effect  on  total  fuel 
consumed  than  the  takeoff  fuel,  with  a  significance  of 
.854  versus  a  negligible  .007  for  the  takeoff  fuel  factor. 

Therefore,  the  null  hypotheses 


V 

“l = 

°2  - 

a3 

(32) 

H  : 
o 

%  - 

02  - 

•  •  •  = 

(33) 

V 

(a8) 

11  = 

(aB) 12  =  • . .  =  (aB) i j 

(34) 

i=l, 2, 3 
j~l ,2,.. .,7 


cannot  be  rejected  at  the  significance  level  of  .95;  hypo¬ 
thesis  (32)  becomes  significant  at  a  level  of  .85,  but  the 
others  are  insignificant.  This  confirms  that  cargo  load 
has  a  much  greater  effect  than  takeoff  fuel  on  total  fuel 
consumed. 

Part  Three 

A  C-5A  mission  from  Rhein  Main  AB,  West  Germany, 
was  selected  to  study  the  effect  of  the  parameters  NSTEP 
and  ITMAX  on  the  model  output.  This  was  done  to  determine 
a  value  for  each  of  these  parameters  that  would  minimize 
computer  execution  time,  since  they  have  the  most  effect  on 
the  total  number  of  computations.  Also,  these  results  were 
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used  with  other  output  during  the  validation  analysis  pre¬ 
sented  in  Chapter  IV.  Table  G-4  contains  the  data  used  in 
this  analysis. 

Other  model  parameters  were  not  changed  from  the 
1  values  recommended  by  Box  (Ref  6) .  These  recommended 

values  were  the  result  of  a  sensitivity  analysis  presented 
*  in  Reference  6,  and  insure  the  most  efficient  application 
of  the  Box  "Complex"  algorithm. 

Part  Four 

Six  typical  aerial  refueling  missions  were  selected 
for  comparison  between  model  output  and  actual  computer 
flight  plan  data.  Wind  factors,  route  of  flight,  tanker 
bases,  and  diversion  bases  were  obtained  from  actual  com¬ 
puter  flight  plans  and  used  as  input  to  the  AR  optimiza¬ 
tion  model.  Since  the  computer  flight  plan  data  were  not 
based  on  optimization  techniques,  the  comparison  shows  the 
AR  optimization  model's  feasibility  and  also  serves  to 
validate  model  output.  Figures  G-3  through  G-8  compare 
the  model  output  for  the  following  airlifter  missions: 

1.  Rhein  Main  AB,  West  Germany  -  Dover  AFB,  DE, 
Figure  G-3  (KC-135A/C-5A) .  A  takeoff  fuel  of  100,000 
pounds  was  used  with  the  model,  versus  225,200  pounds  with 
the  computer  flight  plan.  This  somewhat  distorts  the 
figures  for  comparison,  but  dramatizes  the  effect  of 
excessive  fuel  loads  on  overall  fuel  consumption.  The 
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AIRLIFTER  ROUTE  OF  FLIGHT 


STAGE 
1 
2 

3 

4 

5 

6 

TANKER  DEPARTURE  BASES: 

COORDINATES  BASE 

52.23  -.73  Mildenhall,  England 

47.43  68.36  Loring  AFB ,  ME 

AIRLIFTER  DIVERSION  BASES: 


COORDINATES 

BASE 

50.25 

-8.16 

Rhein 

Main 

,  West  Germany 

39.11 

75.46 

Dover 

AFB, 

DE 

53.31 

60.55 

Goose 

Bay, 

Labrador 

Fig.  G-2.  Input  Data  for  Part  Three 


COORDINATES 

WIND  FACTOR 

DISTANCE 

39.11 

75.46 

-8 

1005.44 

54.81 

66.75 

-32 

764.30 

60.00 

45.00 

-9 

732.29 

62.00 

20.00 

-36 

783.38 

54.11 

.13 

-18 

382.56 

50.25 

-8.16 

0 

0.00 

147 


TABLE  G-4 


NSTEP  AND  ITMAX  EFFECT  ON  OUTPUT 


ITMAX  TOTAL  FUEL  CONSUMED 


NSTEP 


OFFLOAD 


231140 

105684 

228653 

103817 

228748 

107073 

229246 

111154 

229246 

111154 

229246 

111154 

229246 

111154 

229246 

111154 

computer  flight  plan  shows  a  refueling  point  only  382 
nautical  miles  from  that  selected  by  the  model.  Computer 
flight  plan  data  results  in  a  significantly  higher  fuel 
state  overhead  destination  than  required,  also  adding  to 
airlifter  fuel  consumption. 

2.  Tinker  AFB,  OK  -  Kadena  AB,  Japan,  Figure  G-4 
(KC-135A/C-5A) .  Again,  the  computer  flight  plan  calls 

for  an  onload  much  greater  than  that  required,  and  a  sig- 

/ 

nificant  portion  of  this  extra  fuel  is  burned  just  by 
carrying  it.  Fuel  state  values  are  in  gbod  agreement. 

The  AR  point  chosen  was  exactly  the  same  as  that  used  in 
the  computer  flight  plan.  Again,  fuel  overhead  destina¬ 
tion  for  the  computer  flight  plan  is  much  too  high. 


3.  Dover  AFB  -  Moskow,  Russia,  Figure  G-5  (KC-135A/ 
C-5A) .  The  model  calls  for  one  AR  at  stage  3  as  compared 

to  two  ARs  at  stages  4  and  8  for  the  computer  flight  plan. 
From  the  model  results,  it  appears  that  the  first  AR  was 
not  necessary.  Again,  too  much  fuel  offload  and  fuel  over¬ 
head  destination  resulted  in  a  higher  fuel  consumed  figure 
for  the  computer  flight  plan. 

4.  Dover  AFB,  DE  -  Cairo,  Egypt,  Figure  G-6 
(KC-135A/C-5A) .  Again,  the  first  AR  shown  on  the  computer 
flight  plan  was  not  necessary.  This  results  in  the  air- 
lifter  carrying  an  excessive  amount  of  fuel  over  a  large 
portion  of  the  flight.  This  results  in  a  higher  fuel  con¬ 
sumption  and  excessive  fuel  overhead  destination.  AR 
points  and  fuel  states  are  in  good  agreement. 

5.  Tinker  AFB,  OK  -  Hickam  AFB,  HI,  Figure  G-7 
(KC-135A/C-5A) .  This  example  was  included  to  illustrate 
the  strong  agreement  between  computer  flight  plan  fuel 
states  and  those  computed  from  the  model.  The  large  dis¬ 
crepancy  (about  6000  pounds)  at  stage  3  was  caused  by  an 
altitude  change  in  the  computer  flight  plan  for  a  planned 
AR  rendezvous  (no  fuel  offload  was  planned,  however) . 

This  resulted  in  slightly  higher  fuel  consumption.  No  AR 
was  required  here,  nor  was  it  called  for  by  the  model. 

6.  McGuire  AFB,  NJ  -  Bahrain,  Bahrain,  Figure  G-8 
(KC-135A/C-141B) .  This  example  shows  reasonable  agreement 


149 


between  AR  points  and  fuel  states.  Again,  the  computer 
flight  plan  call's  for  excessive  fuel  offloads,  resulting 
in  higher  fuel  consumption  and  higher  than  required  fuel 
overhead  destination. 
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Fig.  G-6 .  Dover  AFB  -  Cairo,  Egypt 
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Fig.  G-8.  McGuire  AFB  -  Bahrain 


ARCP  —  The  Aerial  Refueling  Control  Point  is  the  point 

where  the  tanker  orbits  while  waiting  for  the  airlifter, 
usually  100  nautical  miles  downstream  from  the  ARIP. 

The  tanker  departs  the  orbit  at  the  ARCP  only  when  posi¬ 
tive  radio  contact  is  established  with  the  airlifter 
and  the  airlifter  is  within  80  nautical  miles  of  the 
tanker.  The  ARCP  also  defines  the  point  where  the  two 
aircraft  ideally  complete  their  rendezvous  and  begin 
aerial  refueling. 

ARIP  --  The  Aerial  Refueling  Initial  Point  is  the  first 
point  defining  the  aerial  refueling  track.  The  air¬ 
lifter  must  have  positive  radio  contact  with  the  tanker 
before  leaving  this  point.  Also,  at  the  ARIP  the  air¬ 
lifter  begins  a  descent  to  1000  feet  below  AR  altitude 
and  slows  to  rendezvous  airspeed  in  preparation  for  the 
rendezvous. 

"BUDDY"  AERIAL  REFUELING  —  This  mode  of  aerial  refueling 
occurs  when  the  tanker  and  airlifter  takeoff  in  a  cell 
formation,  and  the  tanker  accompanies  the  airlifter 
along  its  route  of  flight.  The  tanker  offloads  its  fuel 
as  required  by  the  airlifter.  The  tanker  can  then 
either  accompany  the  airlifter  for  another  AR,  or  return 
to  a  recovery  base. 
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CATEGORY  I  ROUTE  —  Simply,  this  refers  to  the  overwater 
portion  of  an  airlift  route  segment.  By  definition, 
the  Category  I  portion  of  a  route  is  that  segment  where 
ground  based  radio  aids  to  navigation  are  inadequate 
to  accurately  position  the  aircraft  at  least  once  each 
hour.  The  definition  requires  the  aircraft  to  fly 
directly  over  the  navigation  aid  to  be  "accurately 
posit ioned. " 

END-AR  POINT  —  This  point  defines  the  end  of  the  aerial 
refueling  track;  if  AR  is  not  complete  by  this  point, 
special  permission  from  Air  Traffic  Control  is  required. 
At  this  point,  both  airlifter  and  tanker  depart  the 
track  and  continue  to  their  destinations. 

GREAT  CIRCLE  ROUTE  —  A  great  circle  route  is  the  shortest 
distance  between  two  points  on  the  earth's  surface.  If 
one  were  to  pass  a  plane  through  the  center  of  the  earth 
and  intersect  two  points  on  the  earth's  surface,  the 
line  defined  by  the  intersection  of  the  earth's  surface 
and  the  plane  would  be  the  great  circle  route  between 
those  two  points. 

INDICATED  AIRSPEED  —  This  is  the  airspeed  read  from  the 
aircraft's  airspeed  indicator.  True  airspeed  is  calcu¬ 
lated  from  indicated  by  applying  corrections  for 


installation  error,  air  compressibility  effects,  tempera 
ture  and  altitude. 

KNOTS  —  Nautical  miles  per  hour. 

M-14  —  M-14  is  a  computer  simulation  model  of  the  Military 
Airlift  Command's  airlift  system.  The  model  makes 
extensive  use  of  probability  theory  and  statistical 
sampling  while  simulating  the  effects  of  uncertainty  on 
the  aircraft  network.  M-14  models  all  MAC  possessed 
aircraft,  crews,  maintenance  functions,  normal  airlift 
missions,  and  aerial  refueling,  to  name  just  a  few  of 
its  many  functions.  Reference  17  explains  the  full 
capability  of  M-14  and  its  limitations. 

MULTIPLE  TANKER  AR  —  As  used  in  this  study,  this  term 
refers  to  the  option  of  refueling  from  more  than  one 
tanker  at  each  rendezvous  point.  This  is  contrasted 
with  the  restriction  that  only  one  tanker  can  be 
accepted  at  a  time;  this  is  a  major  limitation  of  pre¬ 
vious  studies  on  aerial  refueling  optimization. 

NAUTICAL  MILE  —  A  nautical  mile  is  defined  as  1/60  of  a 
degree  of  latitude.  It  is  1.15  statute  miles,  or 


NORMAL  RENDEZVOUS  AERIAL  REFUELING  —  This  mode  of  aerial 
refueling  occurs  when  the  airlifter  and  tanker  takeoff 
separately  and  rendezvous  at  a  predetermined  point  for 
refueling.  This  is  by  far  the  most  common  mode.  The 
tanker  returns  to  a  recovery  base  after  AR  is  complete. 

OPERATING  WEIGHT  —  Aircraft  operating  weight  is  defined 
as  the  total  operational  weight  of  the  aircraft  minus 
fuel  and  cargo  load.  It  includes  the  airframe  weight, 
weight  of  all  internal  stores,  oil,  oxygen,  nitrogen 
(C-5A  only) ,  cargo  handling  equipment,  and  aircrew. 

OPTIMUM  ALTITUDE  —  This  is  a  commonly  used  ferm  to  refer 
to  the  airlifter' s  performance  altitude  ceiling.  This 
altitude  is  defined  as  providing  a  300  feet  per  minute 
rate  of  climb  for  the  C-5A,  or  a  400  feet  per  minute  rate 
of  climb  for  the  C-141B.  It  is  not  an  "optimal"  alti¬ 
tude  in  any  sense. 

TEMPERATURE  DEVIATION  —  Standard  atmospheric  temperature 
is  defined  as  15°  Centigrade  at  sea  level.  To  find  the 
standard  temperature  at  any  altitude,  apply  the  standard 
lapse  rate  of  -2°  Centigrade  per  1000  feet.  Temperature 
deviation  from  standard  is  the  difference  between  the 
actual  outside  air  temperature  ar.d  the  temperature  which 
would  exist  on  a  standard  day.  This  figure  is  used  in 
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lieu  of  the  actual  temperature  during  all  mission  plan¬ 
ning  . 

TRUE  AIRSPEED  —  Intuitively,  true  airspeed  is  the  speed 
of  the  aircraft  through  (or  relative  to)  the  air  mass. 

To  obtain  true  airspeed,  multiply  aircraft  indicated 
Mach  number  by  the  local  speed  of  sound.  Ground  speed 
is  calculated  from  true  airspeed  by  adding  the  wind 
factor. 

WIND  FACTOR  —  Wind  factor  is  defined  as  the  average  effec¬ 
tive  wind  an  aircraft  experiences,  based  on  wind  direc¬ 
tion  and  aircraft  heading,  over  a  route  or  route  segment . 
For  example,  a  wind  factor  of  50  indicates  the  aircraft 
experiences  an  average  tailwind  of  50  knots  over  the 
route  segment.  A  negative  wind  factor  indicates  a 
headwind.  The  wind  factor  is  added  algebraically  to  the 
true  airspeed  to  obtain  ground  speed. 
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With  the  advent  of  strategic  airlift  in-flight  refueling, 
the  role  of  strategic  mobility  in  supporting  our  national  policy 
has  taken  on  new  dimensions.  Airlift  of  important  military 
equipment  to  any  part  of  the  world  is  no  longer  dependent  on 
enroute  bases  for  refueling.  Also,  our  response  time  to  any 
contingency  far  from  home  can  now  be  measured  in  terms  of  hours 
instead  of  days.  This  increased  demand  on  our  limited  aerial 
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refueling  resources,  plus  the  high  cost  and  doubtful  availability 
of  jet  fuel,  require  that  any  deployment  using  aerial  refueling 
be  done  as  efficiently  as  possible. 

The  main  objective  of  this  study  is  to  define  and  develop  an 
operationally  usable  optimization  model  that  will,  given  airlifter 

I  route  of  flight,  takeoff  fuel  and  gross  weight,  enroute.  wind  fac¬ 
tors,  tanker  departure  bases,  and  available  enroute  diversion 
^  bases,  determine  the  aerial  refueling  points,  fuel  offload,  and 
tanker  departure  bases  to  minimize  total  fuel  consumed  by  both  air 
lifter  and  tanker.  The  key  to  this  study  is  the  requirement  that 
the  model  be  operationally  usable.  Therefore,  the  model  is 
designed  to  allow  input  of  any  mission  scenario  using  geographic 
coordinates . 

The  var&njs  subobjectives  presented  in  this  paper  are 
organized  wi'tn  a  view  towards  validating  and  verifying  the  model's 
output.  However,  some  questions  are  also  answered  regarding  the 
nature  of  the  aerial  refueling  optimization  problem.  An  analysis 
of  the  experimental  design  shows  that  it  is  extremely  difficult 
if  not  impossible  to  develop  general  "rules  of  thumb"  to  minimize 
total  fuel  consumed,  primarily  because  of  the  complex  inter¬ 
action  between  the  control  variables  within  the  aerial  refueling 
system.  It  is  found  that  cargo  load,  and  hence  airlifter  gross 
weight,  has  more  effect  on  total  fuel  consumed  than  takeoff  fuel. 
Also,  it  is  found  that  the  objective  of  minimizing  total  fuel  con¬ 
sumed  is  not  always  consistent  with  that  of  minimizing  the  use  of 
tankers. 

The  model  uses  the  dynamic  programming  technique  during  the 
solution  process.  In  addition,  the  "Complex"  Method  of  Box  is 
incorporated  to  optimize,  or  minimize,  stage  return  functions. 

It  is  anticipated  that  the  use  of  this  model  to  plan  opera¬ 
tional  aerial  refueling  missions  will  result  in  significant  fuel 
savings  to  the  Air  Force. 
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